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The  program  consists  of  two  major  task  areas  comprising  an  interactive  numerical/experimental 
approach.  Task  I  involves  the  numerical  computation  of  the  key  SAW  properties  for  doubly  rotated  quartz 
substrates  for  the  purpose  of  locating  promising  angular  ranges  with  properties  superior  to  the  singly 
rotated  cuts  now  in  existence.  More  detailed  calculations  follow  to  refine  the  angular  coordinates  in  order  to 
specify  cuts  for  experimental  verification  in  Task  II.  In  Task  II,  sets  of  substrates  with  promising  orientations 
identified  in  T ask  I  will  be  prepared  and  SAW  device  patterns  will  be  fabricated  for  evaluation  of  the  key  SAW 
properties.  The  experimental  results  of  this  task  will  be  correlated  with  the  theoretical  predictions  and  an 
iterative  process  developed  for  refinement  of  both  theoretical  and  experimental  parameters.  As  the  program 
proceeds,  working  SAW  device  models  are  delivered  as  a  demonstration  of  progress  and  an  indication  of  the 
future  potential  of  the  doubly  rotated  cuts.  This  report  contains  the  results  and  summary  of  Task  I,  Task  II 
and  the  results  of  an  iteration.  The  discovery  of  variation  of  power  flow  angle  as  a  function  of  temperature 
and  methods  for  compensation  of  the  temperature  coefficient  of  frequency  as  the  power  flow  angle  changes 
is  discussed. 
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SECTION  1 
INTRODUCTION 


1.  PROGRAM  OBJECTIVE 

The  objective  of  this  program  is  the  exploratory  development  of  doubly  rotated  cuts  of  quartz 
possessing  superior  Surface  Acoustic  Wave  (SAW)  properties  for  applications  involving  environmentally 
hardened  devices.  The  key  properties  examined  and  optimized  both  theoretically  and  experimentally  are: 
first,  second  and  third  order  Temperature  Coefficients  of  Delay  (TCD),  piezoelectric  coupling  factor,  power 
flow  angle,  Bulk  Acoustic  Wave  (BAW)  inverse  velocity  surfaces,  leaky  waves,  and  sensitivities  of  the  above 
quantities  to  misorientations  and  manufacturing  tolerances. 

2.  PROGRAM  SCOPE 

The  program  consists  of  two  major  task  areas  comprising  an  interactive  numerical/experimental 
approach.  Task  I  involves  the  numerical  computation  of  the  key  SAW  properties  for  doubly  rotated  quartz 
substrates  for  the  purpose  of  locating  promising  angular  ranges  with  properties  superior  to  the  singly 
rotated  cuts  now  in  existence.  More  detailed  calculations  follow  to  refine  the  angular  coordinates  in  order  to 
specify  cuts  for  experimental  verification  in  Task  II.  In  Task  II,  sets  of  substrates  with  promising  orientations 
identified  in  Task  I  are  prepared  and  SAW  device  patterns  fabricated  for  evaluation  of  the  key  SAW 
properties.  The  experimental  results  of  this  task  are  correlated  with  the  theoretical  predictions  and  an 
iterative  process  develops  for  refinement  of  both  theoretical  and  experimental  parameters.  As  the  program 
proceeds,  working  SAW  device  models  are  delivered  as  a  demonstration  of  progress  and  an  indication  of  the 
future  potential  of  the  doubly  rotated  cuts. 

3.  TECHNICAL  APPROACH  SUMMARY  OF  TASK  I 

To  accurately  characterize  the  properties  of  doubly  rotated  quartz,  three  basic  capabilities  are 
essential: 

a.  Theoretical  approach  and  associated  computer  software  which  will  accurately  and  quickly 
locate  promising  zero  TCD  cut  angles  and  characterize  the  other  key  SAW  parameters 

b.  Source  of  rotated  quartz  substrates  of  superior  quality  which  can  be  quickly  fabricated  and 
the  angular  orientation  determined  with  a  high  degree  of  precision 

c.  Required  fabrication  facilities  and  measurement  tools  to  accurately  determine  the  key  SAW 
device  parameters. 
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In  the  theoretical  area,  this  program  has  characterized  two  basic  theoretical  approaches  for  the 
identification  of  zero  TCOs  on  rotated  cuts  of  quartz.  For  this  study,  two  computer  programs  available  at 
Motorola  are  used.  The  first  program  calculates  the  first,  second  and  third  order  TCDs  of  rotated  cuts  using  a 
finite  difference  method.1  This  technique  is  simple,  well  established,  and  has  been  used  for  analytically 
determining  the  temperature  coefficient  curves  for  singly  and  doubly  rotated  cuts  of  quartz.  To  more 
accurately  define  the  temperature  coefficient  properties,  a  second  program  which  encompasses  lattice 
skewing  effects  is  used.  This  more  complete  theoretical  approach  is  based  on  the  work  of  Sinha  and 
Tiersten2,  and  its  utility  has  been  verified. 

The  final  theoretical  work  is  the  characterization  of  the  other  key  parameters  with  standard  SAW 
programs  used  routinely  for  material  characterization  and  device  development. 

Accurately  oriented  quartz  bars,  supplied  by  Motorola,  Carlisle,  are  cut  at  Motorola,  lapped  and 
finely  polished  at  Motorola  and  Crystal  Technology  Inc.  A  mechanical  polishing  procedure  is  used.  During 
this  program,  several  substrates  from  a  single  bar  with  incremental  angular  deviation  about  a  promising 
angular  position  are  fabricated.  By  carefui  organization  of  the  angle  selection  and  cut  procedures,  a 
substantial  savings  in  time  and  money  is  achieved. 

The  angular  orientation  of  the  doubly  rotated  substrates  are  defined  to  an  accuracy  of  within  ±5 
minutes  using  X-ray  diffractometry.  Equipment  used  includes  Laue  pattern  X-ray  equipment,  X-ray 
diffractometers,  and  precision  wafer  cut  and  polishing  equipment. 

4.  TECHNICAL  APPROACH  SUMMARY  OF  TASK  II 

During  this  period  the  first  iteration  theoretical  calculations  performed  to  characterize  doubly  rotated 
cuts  of  quartz  were  completed.  Theoretically  temperature-stable  cuts  with  zeroTCF1"  andTCF12'  as  small  as 
-1 .0  x  10  ®  were  located.  This  represents  a  better  than  three-fold  improvement  over  the  ST  cut.  Experimental 
measurements  of  the  TCF's  have  been  performed  on  some  doubly  rotated  cuts.  Zero  TCF("  SAW  devices  for 
which  the  second  order  temperature  term  is  the  dominant  term  have  been  fabricated.  Measured  values  of 
TCF'21  of  these  devices  as  low  as  -1.5  x  10"8  have  been  obtained.  The  agreement  between  the  experimental 
and  the  calculated  results  was  excellent.  Measured  first  order  TCF's  and  those  calculated  by  Sinha  and 
Tiersten's  perturbation  program  were  found  to  differ  by  less  than  3  ppm/C°,  and  second  order  TCF’s 
calculated  by  the  finite  difference  method  were  found  to  differ  by  less  than  0.005  ppm/C02.  Both  the 
measured  and  calculated  third  order  TCF's  were  found  to  be  too  small  to  be  a  significant  factor  in  device 
performance. 

"‘Numerical  Computation  of  Acoustic  Surface  Waves  in  Layered  Piezoelectric  Media-Computer  Program 
Descriptions”,  William  Jones,  William  Smith,  Donald  Perry,  Final  Report  F 1 9628-70- C-0027 ,  prepared  for  Air 
Force  Cambridge  Research  Laboratories  by  Hughes  Aircraft  Company. 

2"0n  t  he  Temperature  Dependence  of  the  Velocity  of  Surface  Waves  in  Quartz”,  B.K.  Sinha  and  H.F.  Tiersten, 
Proceedings  of  the  32nd  Annual  Symposium  of  Frequency  Control,  1978,  pp.  150-153. 
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A  complete  SAW  test  area  and  optical  laboratory  form  the  basis  for  the  experimental  evaluation  of  the 
key  SAW  parameters  of  the  doubly  rotated  quartz  delay  lines,  oscillators  and  resonators.  The  equipment  is 
set  up  for  rapid  display,  measurement  and  recording  of  propagation  directions,  TCF's,  velocities,  beam 
steering  angles  and  diffraction. 

The  excellent  agreement  between  the  experimental  and  theoretical  results,  and  the  success  with 
which  low  TCF  cuts  have  been  located,  confirmed  the  utility  and  accuracy  of  the  techniques  used  in  the 
program.  The  second  theoretical  and  experimental  iterations  promise  to  yield  orientations  and  devices  with 
even  greater  temperature  stability  than  those  already  obtained. 
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SECTION  II 

TECHNICAL  DISCUSSION  OF  TASK  I 


1.  INTRODUCTION 

Quartz  is  the  most  commonly  used  substrate  for  fabricating  Surface  Acoustic  Wave  (SAW)  devices. 

In  SAW  narrowband  filter,  oscillator,  and  resonator  applications,  the  temperature  stability  of  the  device  is 
an  important  design  parameter.  Currently,  almost  all  SAW  devices  fabricated  on  quartz  use  the  ST-Cut, 
which  exhibits  a  parabolic  frequency  dependence  in  temperature.  For  many  applications,  the  temperature 
dependence  of  devices  fabricated  on  ST  quartz  is  too  large.  Thus  it  is  desirable  to  find  crystal  cuts  with 
superior  temperature  performance.  Of  course,  many  other  design  parameters  must  be  considered  when 
choosing  a  crystal  cut.  Some  of  the  more  important  ones  are  the  piezoelectric  coupling  coefficient,  acoustic 
losses,  dependence  of  device  performance  on  cut  misorientation,  excitation  of  bulk  modes,  and  beam 
steering  angle.  These  parameters  can  all  be  determined  for  a  given  cut. 

The  objectives  of  this  program  are  to  find  crystal  cuts  which  exhibit  zero  temperature  coefficients  of 
delay  so  that  there  will  be  no  frequency-temperature  dependence  observed  in  temperature  stable 
oscillators,  resonators  and  filters.  Computer  models  and  experimental  verification  were  used  to  investigate 
the  temperature  dependence  of  different  cuts  of  crystal  for  SAW  devices. 

Defining  r  as  the  delay  time  for  an  acoustic  wave  to  propagate  between  two  points  on  the  surface  of 
the  crystal,  we  wish  to  find  orientations  for  which  r  is  constant  in  temperature,  or  more  formally  (1/r) 

(dr/dT)  |  I  ;  degrees  ('  —  0. 

If  F  is  the  frequency  of  a  SAW  resonator,  then  the  1st  order  derivative  is:  (1/F)  (dF/dT)  =  -(1/r) 
(dr/dT).  It  is  desirable  to  have  the  higher  order  derivatives  to  be  as  close  to  zero  as  possible,  or  (1/2r) 

(d2r/dT2)  |  T  =  2'  degrees  C  =  0,  6tC. 

Letting  9  be  the  length  between  two  points,  (r)  is  simply  given  by  r  =  (2/Vs),  (1/r)(dr/dT)  = 
(1/8)(d2/dT)  -  (1/V,)(dV,/dT). 

Two  computer  programs  for  calculating  r  as  a  function  of  the  stress,  the  dielectric  and  piezoelectric 
constants  of  a  substrate  material  are  available  for  this  study.  Furthermore,  the  temperature  variation  of 
those  constants  for  quartz  are  available,  allowing  one  to  calculate  (1/V.,(dVs/dT)  by  a  finite  difference 
method  (note  that  higher  order  terms  can  be  calculated  the  same  way).  If  the  crystal  expansion  coefficients 
are  known,  it  is  then  a  simple  matter  to  calculate  (1/r)(dr/dT)  for  any  particular  orientation.  Perturbation 


programs  developed  by  Sinha  and  Tiersten  to  calculate  (1/T)(dr/dT)  are  also  available.  Combined  with  a 
search  method,  the  cuts  for  which  first  order  temperature  coefficient  of  frequency  vanishes  can  be  found. 

a.  Calculation  of  Temperature  Coefficients 

It  has  been  shown  that  determining  the  temperature  dependence  of  r  (time  delay)  is  equivalent 
to  determining  the  temperature  dependence  of  F  (frequency)  via  the  relation  F  °°  1/t. 

Since  the  experimental  data  is  derived  from  frequency  measurements,  the  frequency 
characteristics  of  the  devices  is  first  determined.  The  relation  between  the  temperature  coefficient  of 
frequency  (TCF)  and  temperature  coefficient  of  delay  (TCD)  are  related  as  follows  (See  Appendix  E): 


a  f 

«lf' 

=  -«^  +  («r' 

a 131 

=  -a'3,+  2a' 

where  an"’  is  the  ith  order  TCF,  a  T is  the  ith  order  TCD. 

Using  the  relations  above,  one  can  always  relate  one  set  of  the  temperature  coefficients  to  the 


The  various  procedures  outlined  in  the  following  sections  will  yield  the  quantities  a,"',  i  =  1,2, 
3,  with  V,(T)  3  V,(T„)(1  +  «,"’  dT  +  a,'2’  dT2  +  a,'3'  dT3)  and  dT  =  T  -  T„. 

The  problem  of  interest  is  not  only,  however,  in  finding  the  a,(”’s  but  in  finding  the  delay  time  r 
and  the  frequency  dependence  F  of  a  device.  The  frequency  dependence  F(T)  =  F(T0)(1  +  aF(”  dT  +  aF12’  dT2 + 
ara'  dT3)  is  a  function  of  not  only  V,  but  of  n(T0)  =  2(T0)(1  +  a,"’  dT  +  a,'2’  dT2  +  a,'31  dT3),  the  spacing  between 
reflectors  in  a  resonator  (or  the  length  of  the  delay  line  in  an  oscillator) 

From  the  above  relations  and  (1  +  X)-’  a  1  -  X  +  X2  -  X3  if  X  «1,  then: 

F  =  V,  2  -  f  V,  (1  +  «v‘"  dT  +  W21  dT2  +  av'3’  dT3)}/2„  (1  +  a,"3  dT  +  a,'2'  dT2  +  a,'3'  dT3)  (4) 

as  F(T„)  (  1  +  av"’  dT  +  av12’  dT2  +  av'3’  dT3  {5] 

-  a,"’  dT  -  a,"’  av"'  dT2  -  a,"'  av'2’  dT3 

-  a,'21  dT2  -  a,'2’  a,"'  dT3  -  a,'3’  dT3 

+  (a,"’)2  dT2  +  (a,"')2  av"’  dT3 

+  2a,"’  a,'2’  dT3  -  (a,"’  dT)3}  (6 


(7) 

(8) 
(9) 


a»  —  0\  —  os  a\  +  [on  f  -  ai 

-  —  »'»  ^o>  i-)  i  n»  /  n»v2 

a t  —  a4  —  ai  a*  +  a*  (ai  ) 

— )  +  4Qi  on  —  a,  ai  —  ai 

giving  the  frequency  dependence  directly.  The  calculation  of  the  temperature  coefficients  of  velocity,  is 

achieved  by  calculating  the  phase  velocity  with  the  Finite  Difference  Technique  or  other  techniques 
discussed  below  for  a  variety  of  temperatures.  As  the  velocity  is  a  function  of  temperature,  a  linear 
regression  program  is  used  in  the  finite  difference  technique  to  curve  fit  the  data  to  a  third  order 
polynominal.  The  constants  a,l2>  and  a*0'  are  thus  obtained  by  optimum  curve  fitting  of  the  data  points 
to  the  polynominal.  The  temperature  coefficients  of  length,  are  found  in  standard  references.  It  should  be 

noted  that  the  a,’s  also  depend  on  the  direction  of  propagation,  and  must  be  calculated  for  each  direction  of 
propagation  considered  by  a  simple  geometrical  transformation  analogous  to  that  used  to  rotate  all  of  the 
other  physical  constants. 

b.  Rayleigh  Wave 

We  will  briefly  discuss  the  Rayleigh  wave  solutions  in  their  general  form  in  this  section.  The 
coordinate  system  is  defined  with  the  Z  or  1  axis  being  the  direction  of  propagation  and  the  -Y  or  3  axis 
normal  to  the  crystal  surface.  Vs  will  denote  the  Rayleigh  wave  velocity,  /SR = a>/V„  the  wave  number,  u,  (i = 1 , 
2, 3),  the  particle  displacements  along  the  1, 2,  or  3  axis,  u4  =  <t>,  the  electric  potential.  Also,  T„  denotes  the 
stress  tensor. 

Sij  =  (1/2)(du,/dx,  +  dUj/dx,),  the  strain  tensor,  cijkl,  the  elastic  tensor,  p  the  density  of  the 
substrate,  D,  the  electric  displacement,  E,  the  electric  field,  the  dielectric  constant,  eijk  the  piezoelectric 
constant.  We  also  use  a  for  the  difference,  where,  for  example, 

Ap  =  p(T)  -  p(To)  =  p(T0)(a  p,0dT  +  a  p,2idT2  +  a  pi3)  dT3),  and  let  P  =  power/unit  width  in  the 
x  direction. 

We  assume  relations  such  as 

T  =  — e  E  +  cE:S  or  S  *  e-E  +  SE:T 
0  =  «SE  +  e:S  0  =  «TE  +  d:T 
pd2  u,/at2  -  Cijn  a2uk/ax,  ax,  -  ektj  atyax,  axk  =  0 
e,u  a:uk/ax,  ax,  -  «,k  a2  <t>i  ax,  axk  =  0 


(10) 

(ID 
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The  general  solution  for  a  wave  traveling  on  the  surface  can  be  written 


u,  =  [  £  C™  a,w  exp(-(i<o/V.)  ii)  b‘r*  x()  ]  exp(iwt)  for  j  =  1  to  4  (12) 

m  1*1 

4 

For  the  Rayleigh  wave,  this  reduces  to  U,  =  £  cm  «/->  exp  H/Sub^y)  exp(i/3Rz  -  i«t)  (13) 

m- 1 


where  the  coefficient  cm,  ^""’(weighting  factors),  /3R(wave  number),  and  b<m|(complex  decay  constants)  are 
to  be  calculated  by  the  standard  iterative  procedure  on  a  computer. 

The  fields  are  calculated  as  E  =  Vu4  with  E,  =  0,  Ey  =  (-i0Rb<m,)u4  and  E*  =  (i0R)u4. 

c.  Methods  for  Calculating  the  Temperature  Dependence  of  the  Rayleigh  Velocity 

In  the  following  sections,  four  methods  for  calculating  the  temperature  characteristics  of  the 
Rayleigh  wave  are  discussed. 


(1)  Finite  Difference  Technique  of  Calculating  Temperature  Dependence  of  the  Rayleigh 

Velocity 


Before  TC0("  =  (1/r)(dr/dT)  can  be  calculated,  one  must  first  calculate  the  dependence 
of  the  Rayleigh  wave  velocity  on  temperature.  The  most  straightforward  method  for  doing  this  is  the  finite 
difference  method.  The  Rayleigh  wave  velocities  are  calculated  for  different  temperatures,  yielding  the 

values  V,(T,),  i  ==  1, 2 n.  This  is  done  by  first  calculating  the  fundamental  constants  at  the  temperature  T,  of 

interest.  The  fundamental  constants  are  then  rotated  into  the  coordinate  system  of  interest.  An  iterative 
procedure  is  used  to  calculate  a  velocity  V,  for  which  Christoffel’s  equation  and  the  boundary  conditions 
are  satisfied  simultaneously  (see  Appendix  D).  Simple  finite  difference  techniques  can  be  used  to  calculate 
(1  /V,)(dV,/dT),  (1  /2V,)(d2V,/dTJ),  etc.  For  example,  after  calculating  V,(T.)  for  T,  =  T„,  L  =  T0  +  AT,  L  =  To  - 
AT,  we  can  use  (1  / V,)(dV,/dT)  at  Tn  ss  (1  /V0(To))(V,(T0  +  AT)  -  V,(T0  -  aT))/2aT.  Alternately,  standard  linear 
regression  of  polynomials  may  be  employed  to  yield  those  coefficients.  Another  approach  consists  of 
calculating  directly  the  frequency-temperature  characteristics  of  the  orientation  for  several  widely  differing 
temperatures.  A  measure  of  the  temperature  stability  is  then  used.  Letting  F(T)  be  the  frequency  of  the 
device,  the  measure  of  deviation  is  calculated  as 


RMS  frequency  deviation  = 


[(X(F(Ti)-F(T0))2)/n-  F  (T„)2]l/2  with 

i  *  I 

F(T)  =  V.  (T)/20(1  +  aL"’(AT)  +  aL'2,(AT)2  +  aL°'(AT)3), 


(14) 


although  a,"',  on'21,  and  of"1  are  calculated  for  purposes  of  comparison. 
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In  much  of  the  earlier  work,  the  SAW  velocities  were  calculated  at  three  temperature 
points,  T,  =  -50°C,  T: = 25°C,  and  T,  =  100°C,  to  save  on  computer  time.  The  following  formulas  were  used  for 
the  calculations  of  first  and  second  order  temperature  coefficients  of  velocity. 


=  [V(T ,)  -  V(Tj)]/[V(T-’)  (T,  -  T,)]  (15) 

=  [V(100)  -  V (—50)]  / (V (25)  •  150] 

oj  =  [V(T ■)  +  V(T,)  -  2  •  V(L)]/[V(T2)  •  2  (T3  -  L)(T2  -  T,)J  (16) 

=  [V(100  +  V(— 50)  -  2  •  V(25)]/[V(25)  ■  2  (75)(75)] 

In  the  more  recent  works,  six  temperature  points  were  calculated  for  each  orientation  to  provide  data  for  the 
linear  regression  analysis.  The  six  temperature  points  were  100°C,  85°C,  75°C,  -25°C,  -10°C,  and  ~50°C.  The 
results  were  used  to  calculate  TCF"',  TCF'21  and  TCF"'. 

(2)  Perturbation  Technique  of  Calculating  Temperature  Dependence  of  the  Rayleigh  Velocity 

Pertubation  theory  may  be  applied  to  the  problem  of  calculating  the  first,  second  and 
third  order  dependence  of  the  Rayleigh  velocity  V,  of  a  piezoelectric  substrate,  once  the  solution  to  the 
Rayleigh  wave  propagation  at  a  temperature  reference  T0  is  known. 

Pertubation  techniques  allow  calculations  of  small  changes  in  the  solutions  to  a  problem 
caused  by  small  changes  in  the  physical  parameters  of  the  problem,  once  the  solution  to  the  unperturbed 
problem  is  known.  One  can  apply  perturbation  techniques  either  to  boundary  perturbation  such  as  mass 
loading  on  the  surface  or  to  volume  perturbations  such  as  adding  a  thin  conducting  layer  to  the  interior  of 
piezoelectric  substrates.  As  the  boundary  conditions  are  unaffected  by  changes  in  temperature  while 
material  constants  such  as  c„  are  temperature  dependent,  one  can  apply  the  volume  perturbation  theory  to 
the  problem. 

The  general  approach  to  the  problem  of  determining  the  temperature  dependence  of  V, 
will  be  as  follows.  First,  the  Rayleigh  wave  propagation  problem  will  be  solved  in  the  standard  way  in  its 
entirety  at  room  temperature,  T0.  Given  the  solution  of  the  problem  at  T0  and  the  dependence  of  the  physical 
constants  (such  as  c„)  on  temperature  at  T0,  one  will  apply  the  volume  perturbation  formula,  calculating  the 
temperature  dependence  of  V„  The  dependence  of  V,  on  T  is  then  used  to  calculate  the  frequency 
characteristics  of  the  actual  device  given  the  thermal  expansion  coefficients  as  a  function  of  temperature.  At 
this  point,  the  frequency  temperature  dependence  of  the  substrate  as  a  function  of  the  crystal  cut  and 
direction  can  be  thoroughly  explored.  See  Appendix  A  for  a  complete  discussion  of  this  method. 
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(3)  Approach  of  Sinha  and  Tiersten 


The  primary  difficulty  with  the  perturbation  technique  is  that  it  does  not  take  into 
account  the  change  of  coordinate  systems  induced  by  the  thermal  expansion  in  the  material.  This  comes 
about  because  the  material  distorts  as  temperature  changes.  Thus,  the  set  of  axes  to  which  the 
fundamental  elastic  constants  refer,  which  is  fixed  to  the  crystal,  is  no  longer  equivalent  to  the  axes  used  to 
calculate  V,.  This  problem  is  elegantly  solved  by  Sinha  and  Tiersten1*2.  The  first  simplification  which  occurs 
is  that  the  density  of  the  material  remains  constant  with  temperature.  Furthermore,  the  aF("  simply  become 
oi =  a,1",  as  -f-a  length  in  this  coordinate  system  does  not  change.  The  only  difficulty  is  that  the  elastic 
constants  previously  used  no  longer  refer  to  the  proper  coordinate  system  and  must  be  rederived  from  the 
original  experimental  data.  This  procedure  has  already  been  carried  out  for  the  first  order  temperature 
derivation  of  quartz2  but  has  not  deen  done  for  the  second  and  third  order  coefficients.  Nonetheless,  the 
procedure  yields  more  accurate  results  for  the  first  order  dependence.  See  Appendix  B  for  a  mathematical 
description  of  the  salient  features  of  this  technique. 

(4)  Differentiation  Method 

A  method  for  determining  the  theoretical  temperature  dependence  of  Rayleigh  Surface 
Waves  consists  of  formally  differentiating  the  wave  equation  and  boundary  conditions.  The  boundary 
conditions  and  wave  equation  must  be  true  at  alt  temperatures,  placing  constraints  on  how  the  parameters 
of  the  wave  equation  may  vary.  In  this  technique,  the  derivatives  of  these  equations  with  respect  to 
temperature  are  set  to  zero  and  solved  for  the  velocity-temperature  dependence.  This  method  follows  the 
methods  used  by  Bechmann,  Ballato,  and  Lukaszek  to  compute  the  temperature  dependence  of  the 
fundamental  elastic  constants  from  frequency  data,  except  that  the  simplifying  assumptions  of  assuming 
bulk  wave  solutions  cannot  be  made.  This  method  was  later  used  by  Hauden’  to  search  for  temperature 
stable  cuts  of  quartz.  See  Appendix  C  for  a  discussion  of  this  method. 

(5)  Summary  of  Approaches 

Of  all  the  techniques  presented,  the  finite  difference  technique  satisfies  all  of  the  basic 
requirements  for  calculating  temperature  coefficients  of  delay.  Arbitrary  crystal  structures  may  be 


'“On  the  Temperature  Dependence  of  the  Velocity  of  Surface  Waves  in  Quartz,"  Sinha  and  Tiersten.  1978 
Ultrasonics  Symposium  Proceedings,  pp.  662-665. 

^‘Temperature  Dependence  of  the  Fundamental  Elastic  Constants  of  Quartz,”  Sinha  and  Tiersten, 
Proceedings  of  the  32nd  Annual  Symposium  on  Frequency  Control,  1978,  pp.  150-153. 

’"Higher  Order  Temperature  Coefficients  of  Quartz  SAW  Oscillators,"  D.  Hauden,  M.  Michael,  J.  J.  Gagnepain, 
Frequency  Control  Symposium  (1978),  pp.  77-86. 
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investigated  once  the  density,  piezoelectric,  elastic,  and  dielectric  constants  and  their  temperature 
variation  are  known.  When  double  precision  on  the  computer  is  used  and  when  the  velocity  over  a  large 
temperature  range  is  calculated,  it  becomes  a  very  accurate  numerical  approach.  The  computational 
efficiency  of  the  finite  difference  method  is  not  as  high  as  some  of  the  other  techniques,  owing  to  the  need  for 
repeated  calculation  of  the  surface  wave  velocity. 

A  comparison  between  the  experimental  results  of  Schulz1  and  Motorola's  finite 
difference  program  adapted  from  Jones  et  al.2  as  shown  in  Figure  1.  The  results  for  this  cut  and  many  other 
orientations  studied  have  been  found  to  be  very  good.  The  slight  shift  of  the  analytical  curve  versus 
experimental  curve  can  be  explained  by  a  crystal  misorientation. 


Figure  1.  Finite  Difference  Method  Versus  Experiment  ST-Cut  (Quartz) 


'“Surface  Acoustic  Wave  Delay  Lines  with  Small  Temperature  Coefficient,”  Schulz,  Manfred  B.,  Proc.  IEEE, 
Sept.  1970,  pp.  1361-1362. 

2“Numerical  Computation  of  Acoustic  Surface  Waves  in  Layered  Piezoelectric  Media  -  Computer  Program 
Descriptions",  Jones,  Smith,  and  Perry,  Hughes  Aircraft  Company,  Final  Report,  Air  Force  Cambridge 
Research  Laboratories,  Contract  No.  F19628-70-C-0027. 
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Figure  2  displays  a  phase  delay  versus  temperature  plot  for  the  orientation  (YXwIt) 
56/40/27  in  quartz  obtained  from  Hauden's  graphs2.  Using  the  differentiation  approach,  he  calculated  on-" '  = 
0,  at'2'  =  -0.04  x  10"J  ppm/0C^v>^|,|  =  -22.7  x  10"6  ppm/°C3  for  a  cut  close  to  (YXwIt)  56/40/27.  The  Finite 
Difference  approach  was  used  to  verify  his  first  order  calculation  but  is  in  disagreement  on  the  second  and 
third  order  calculations.  This  can  be  resolved  by  experiment. 


Tiersten’s  Perturbation  program,  while  very  difficult  to  extend  to  higher  orders,  is 
available  for  calculating  the  first  order  TCD. 

Figure  3  shows  the  results  of  Tiersten's  calculation1  and  compares  those  results  with 
the  Finite  Difference  method  as  well  as  experiment.  The  cut  used  for  comparison  is  the  AT-Cut.  The 
improvement  in  accuracy  is  substantial.  It  offers  the  significant  advantages  of  being  quite  cost  efficient,  as 
well  as  being  more  accurate  in  the  first  order  on  off-axis  cuts,  our  primary  area  of  interest.  In  the  task  of 
making  quartz  cuts  along  the  surface  of  zero  TCD,  accurate  values  for  the  first  order  TCD  prove  more  useful 
than  the  less  exact  values  for  the  first,  second,  and  third  order  TCD’s  provided  by  other  methods,  because  it 
is  a  necessary  (but  not  sufficient)  condition  for  zero  TCD  cut  to  get  a  zero  first  order.  Since  it  is  the 
dominating  term,  it  is  important  to  first  get  an  accurate  first  order  zero  TCD  locus.  Then  one  may  search  for 
the  intersection  of  this  locus  with  the  zero  second  and  third  order  terms  with  the  finite  difference  technique. 

'"Temperature  Dependence  of  the  Fundamental  Elastic  Constants  of  Quartz,”  B.  K.  Sinha  and  H.  F.  Tiersten, 
Proceedings  of  the  32nd  Annual  Symposium  on  Frequency  Control,  1978,  pp.  150-153. 

^'Higher  Order  Temperature  Coefficients  of  Quartz  SAW  Oscillators,  ”D.  Hauden,  M.  Michael,  J.  J.  Gagnepain, 
Frequency  Control  Symposium  (1978),  pp.  77-86. 
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9  (DEGREES) 


TEMPERATURE  COEFFICIENTS  OF  ACTUAL 
VELOCITY  ANO  OELAY  FOR  SURFACE  WAVES  ON 
AT-CUT  QUARTZ  AS  A  FUNCTION  OF  PROPAGA¬ 
TION  DIRECTION  RELATIVE  TO  THE  DIAGONAL 
AXIS  AT  2S°C.  THE  DOTTED  CURVE  SHOWS 
THE  AVERAGE  OF  THE  CALCULATED  VALUES  AT 
0°C  AND  50°C  FROM  REF  1.  THE  CIRCLES 
ARE  THE  AVERAGE  OF  THE  EXPERIMENTAL 
VALUES  AT  (PC  ANO  50°C 
GRAPH  OBTAINED  FROM  REF  1 


Figure  3.  Tiersten’s  Method  Versus  Finite  Difference  and  Experiment 


(6)  Investigative  Approach 


A  necessary  but  not  sufficient  condition  to  find  a  temperature  stable  cut  of  quartz  is  that 
ok1"  =  0.  In  practice,  a  sufficient  condition  for  finding  a  zero  temperature  cut  is  that  ap1"  =  of121  =  <*f131  =  0. 
Thus  every  zero  temperature  cut  must  be  on  the  locus  of  angles  which  satisfy  the  condition  ap"’ = 0.  Thus  the 
first  problem  is  to  locate  accurately  such  cuts.  Both  the  finite  difference  approach  and  Tiersten’s  method 
were  used  to  calculate  the  TCF. 


The  approach  used  in  this  program  is  to  identify  the  areas  (angles)  where  TCF(1),  TCFU) 
and  TCF,J)  are  relatively  low  by  using  the  Finite  Difference  Approach  and  then  define  the  exact  orientation 
that  has  zero  TCF"’  in  those  areas  by  the  Sinha  and  Tiersten  approach. 

In  case  the  condition  TCF"’=  TCFI2) = TCF13’ = 0  cannot  be  met,  a  compromise  approach 
would  be  to  find  an  orientation  where  the  first  and  third  order  effects  tend  to  cancel  out  in  the  temperature 
range  of  interest  and  to  find  the  minimum  second  order  effect  orientation  amount  for  those  cuts.  This 
approach  would  provide  an  effective  low  TCF  cut  of  quartz  for  SAW  application. 
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(7)  Analytical  Results  on  Zero  TCF  on  Quartz 


IRE  standard  angle  definitions  (YX  wit)  PHI/THETA/PSI  for  quartz  were  used  throughout 
the  investigation1.  Consider  the  TCFs  to  be  functions  of  these  angles,  which  define  an  angular  volume  0  ^  PHI 
<  30°,  -90°  ^  THETA  ^  90°,  0  PSI  180°, which  spans  the  space  of  possible  cuts  and  propagation 
directions.  The  set  of  points  at  which  TCP"  (PHI,  THETA,  PSI)  =  0  forms  a  surface  in  this  angular  volume. 
Likewise,  the  set  of  points  at  which  TCF'2’  (PHI,  THETA,  PSI)  =0  also  form  surfaces  in  this  angular  volume. 


Our  object  is  to  find  a  point  where  TCF'"  =  TCF'2’ =TCF13’  =  0.  If  the  surface  of  zero  TCF"’ 
intersected  with  the  surface  of  zero  TCF12’,  the  result  would  be  a  line  (or  a  point  if  the  two  surfaces  are  tangent 
to  each  other)  of  angular  points  on  which  TCF1"  =  TCP21 = 0.  The  intersection  of  this  line  with  the  surface  on 
which  TCF15’  =  0  would  yield  a  single  point  at  which  TCP" = TCP2’ = TCF13’ =0.  Neglecting  higher  order  terms, 
we  would  have  found  a  temperature  stable  cut. 

The  calculated  values  of  TCP"  versus  propagation  angles  are  shown  in  Figure  4.  The 
zero  TCP"  is  identified  by  the  areas  where  TCP"  changes  sign. 


Using  the  Finite  Difference  approach  with  the  available  crystal  constants,  the  calculated 
results  show  that  the  zero  TCP"  surfaces  do  not  intersect  with  the  zero  TCP2’  surfaces,  based  on  the 
interpolated  results  of  the  10°  x  10°  x  10°  resolution.  It  is  not  likely  that  a  finer  resolution  will  provide 
contrary  information  because  TCP"  and  TCP2’  are  relatively  slow  varying  functions  as  shown  in  Figure  4. 


'Standards  on  Piezoelectric  Crystals  1949,"  Proc.  IRE  14,  Dec.  1949,  pp.  1378-1395. 
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Calculations  were  performed  on  a  10°  x  10°  x  10°  grid  over  the  angular  ranges  0^  PHI  ^ 
30°,  0  PSI  180°,  and  -90°  THETA  <  90°.  These  angular  ranges,  due  to  the  symmetry  of  quartz,  include 
all  possible  angular  orientations.  These  initial  calculations  defined  the  “angular  volumes”  of  low  TCP 
orientations.  Calculations  were  then  performed  on  a  2.5°  x  2.5°  x  2.5°  grid  near  promising  orientations.  In 
this  way,  the  entire  angular  range  was  explored  and  a  large  computer-based  data  file  built.  Maps  of  first  and 
second  order  TCF’s  are  shown  in  Figure  5.  Zero  first  order  TCF  contours  are  drawn.  In  addition,  contours  of 
the  second  order  TCF  are  drawn  and  shaded  in  areas  where  the  second  order  TCF  is  less  than  0.01  PPM/°C2. 
This  represents  a  substantial  improvement  over  ST  quartz,  for  which  the  second  order  for  TCF  is 
approximately  0.03  PPM/°C2.  Where  data  points  were  missing,  lines  were  connected  by  interpolation. 
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Figure  5.  TCF  Map  of  SAW  Device  on  Quartz 

Despite  the  number  and  density  of  points  at  which  the  first  and  second  order  TCFs  were 
calculated,  wherever  TCF"’  was  found  to  be  less  than  or  equal  to  zero,  we  found  TCF12’  to  be  less  than  0.  The 
reason  for  this  probably  lies  in  the  lack  of  independence  of  the  crystal  constants  themselves.  Using  the 
volume  perturbation  approach  of  Auld1,  one  obtains  equation  A-8  of  Appendix  A. 


"‘Acoustic  Fields  and  Waves  in  Solids,"  Auld,  B.  A.,  Vol.  II,  John  Wiley  &  Sons,  1973,  N.  Y..  p.  297. 
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For  small  changes  in  p  and  cu',  one  might  expect  that  the  TCFs  are  well  correlated  with  the  temperature 
derivatives  of  the  fundamental  constants  for  Rayleigh  waves  in  quartz.  A  list  of  some  of  the  crystal  elastic 


constants’1  temperature  derivatives  is  in  Table  1. 

TABLE  1.  CRYSTAL  ELASTIC  CONSTANTS’  TEMPERATURE  DERIVATIVES 


TC(" 

TC'2> 

TC13’ 

Cm 

-48.5  xIO6 

-107  xIO'9 

-70x10": 

Cu 

-3000 

-3050 

-1260 

Cm 

-550 

-1150 

-750 

Cu 

101 

-48 

-590 

Cu 

-160 

-275 

-250 

C44 

-177 

-216 

-216 

Cm* 

178 

118 

21 

p 

-34.92 

-15.9 

5.3 

an 

13.7 

6.5 

-1.9 

au 

7.5 

2.9 

-1.5 

Figure  6  shows  a  plot  of  these  values.  The  correlation  is  quite  high,  except  for  the  case  of 
Cu  and  suggests  that  the  TCFs  should  be  correlated  also.  Plots  of  TCP"  and  TCP2' versus  angles  shown  in 
Figure  4  reflect  this  correlation. 

Of  course,  if  the  two  TCF  surfaces  do  not  intersect  to  form  a  line  of  zero  TCP"  and  TCP21, 
we  must  use  a  slightly  different  approach  to  finding  temperature  stable  cuts.  The  largest  variation  in 
constants  occurs  for  Cu.  As  we  are  interested  in  the  temperature  range  from  -50°C  to  100°C,  the  maximum 
deviation  of  T  from  To  =  25°C  is  ( T  — T«|  =  75°C.  Calculating  the  changes  in  the  constant  Cu  for  this  maximum 
temperature  change  gives: 


1 

~C^ 

1  A  Cu  1  = 

2.25  x  10- 

1 

Cu 

1  A  Cu"  |  = 

1.72  x  10 

1 

Cu 

(31 

1  A  Cu  |  = 

5.32  x  10" 

(17) 


"‘Higher  Order  Temperature  Coefficients  of  the  Elastic  Stiffnesses  and  Compliances  of  Alpha-Quartz", 
Bechmann,  Ballato,  Lukaszek,  Proc.  IRE,  Aug  1962,  pp.  1812-1922. 
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Thus  the  most  significant  term  is  the  first  order  term.  The  second  order  term  is  still  quite  significant  but  one 
order  of  magnitude  below  the  first,  while  the  third  order  term  is  almost  three  orders  of  magnitude  below  the 
first.  Since  the  delay  time  r  is  dependent  on  the  crystal  constants,  the  most  significant  term  should  be  TCP11. 
The  least  significant  term  should  be  TCP3'.  After  numerous  calculations  of  the  TCFs  for  many  orientations, 
this  appears  to  be  the  case. 

The  investigative  approach  used  has  been  to  first  locate  the  surfaces  of  zero  TCP"  (the 
most  significant  term)  with  the  finite  difference  program.  Near  these  surfaces  of  zero  TCP",  low  values  of 
TCP2'  are  sought,  using  already  calculated  results  of  the  finite  difference  programs.  Where  low  values  of 
TCP2’  have  been  found,  the  perturbation  approach  was  used  to  more  accurately  locate  the  zero  TCP" 
surface,  this  being  the  most  significant  term  in  the  total  temperature  dependence.  TCP3’’s  are  then 
calculated  to  assure  that  their  effect  on  the  total  temperature  dependence  is  small.  To  date,  this  has  always 
been  found  to  be  the  case. 

(8)  Results  of  the  Investigative  Approach 

Table  2  consists  of  a  summary  of  the  results  of  using  the  investigative  approach 
described  in  paragraph  (7)  above.  There  are  three  areas  where  low  TCF cuts  have  been  located.  These  areas 
are  centered  near  (YX  wit)  0/27/138,  (YX  wit)  7/27/135.5,  and  (YX  wit)  15/40/40.  These  orientations  have 
zero  TCP",  calculated  by  the  Sinha  and  Tiersten  approach,  with  TCP2’  and  TCP3’  calculated  with  Finite 
Difference  approach.  These  areas  were  chosen  because  of  zero  TCP",  low  TCP2’,  and  a  low  TCP3’  which  can 
be  mostly  cancelled  out  by  the  first  order  TCF  if  the  propagation  direction  is  slightly  rotated  away  from  the 
zero  TCP"  direction.  Then  the  TCP2’  term  will  dominate  the  performance  characteristics.  The  angular 
resolution  in  these  areas  is  1°  x  1°  x  1°.  The  cuts  potentially  have  one  half  to  one  third  the  temperature 
coefficients  of  ST-Cut  quartz. 
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TABLE  2.  PROPAGATION  CHARACTERISTICS  OF  SELECTED  ORIENTATIONS 


Angles  of  ZTCF'", 

Degrees 

(S  &  T’s  Program) 

TCF(,,/°CJ  (x10-‘) 

Finite  Difference 

Program 

TCF'3'/°C5  (xIO10) 

Finite  Difference 

Program 

Phi 

Theta 

Psi 

6 

26 

136.31 

-1.4 

6 

27 

135.93 

-1.3 

0.67 

6 

28 

135.59 

-1.3 

0.57 

■  1 

26 

135.99 

-1.5 

B 

27 

135.64 

-1.4 

28 

135.27 

-1.3 

0.65 

8 

26 

135.74 

-1.4 

0.65 

8 

27 

135.36 

-1.4 

8 

28 

134,97 

-1.3 

1 

26 

137.78 

-1.2 

0.68 

1 

27 

137.48 

-1.2 

0.65 

1 

28 

137.17 

-1.1 

0.67 

0 

26 

138.07 

-1.2 

0.67 

0 

27 

137.78 

-1.1 

0.68 

0 

28 

137.49 

-1.1 

0.62 

-1 

26 

138.37 

-1.2 

0.60 

-1 

27 

138.09 

-1.2 

0.62 

-1 

28 

137.80 

-1.1 

0.73 

14 

39 

40.195 

-1.0 

0.64 

14 

40 

40.415 

-1.0 

0.66 

14 

41 

40.64 

-1.0 

0.75 

15 

39 

39.79 

-1.0 

0.63 

15 

40 

40 

-1.0 

0.74 

15 

41 

40.23 

-1.0 

0.73 

16 

39 

39.4 

-1.0 

0.68 

16 

40 

39.605 

-1.0 

0.66 

16 

41 

39.825 

-1.1 

0.60 

d.  Piezoelectric  Coupling  Factor 

The  piezoelectric  coupling  factor,  denoted  by  k2,  is  a  measure  of  the  coupling  efficiency  for  an 
interdigitated  electrode. 

We  can  express  k2  in  terms  of  the  interdigital  transducer’s  input  conductance1  as 

G,„M  =  -7T2  6u„  M2Wt.,oe"2  G,P„-,2  (cos  r,  tt)/K2  (cos  t,  tr/2)  (18) 

where  <u„  =  2ir  (2n  -  1)/\s0,  M  is  the  number  of  electrode  pairs,  w  is  the  length  of  the  electrodes,  x  is  the 
spatial  periodicity  of  the  array,  the  width  of  each  electrode  is  17X/2,  is  the  value  of  the  effective 
permittivity  at  zero  velocity,  K  is  the  complete  elliptic  integral  of  the  first  kind,  and  Pm  are  the  Legendre 
polynomials.  G*  is  the  residue  of  1/(  |s  |e.). 

k?  is  defined  by 

k2  =  -2*'°°'  G,  (19) 

Thus  (18)  becomes 

G,„  (<u„)  =  (1/2)  7T2  Oin  M2  W*,100'  k,2  PB-,2  (COS  r,  »r)/K2  (COS  r,  tt/2)  (20) 

ana  we  see  that  k2  is  related  to  the  input  conductance  and  hence  to  the  coupling  efficiency.  A  good 
approximation  to  k\  is  given  by 

k,2  a  2(V„  -  V,)/V,  =  2AV/V  (21) 

where  V,  is  the  shorted  Rayleigh  wave  velocity  and  V„  is  the  free  Rayleigh  wave  velocity. 

It  is  seen  that  the  quantity  aV/V  is  an  important  parameter  as  it  is  a  direct  measure  of  the 
coupling  efficiency  and  conductance  of  an  interdigitated  electrode  pair  on  a  piezoelectric  substrate. 
Furthermore,  aV/V  is  easily  determined  once  the  shorted  and  free  surface  problems  have  been  solved. 
Programs  are  available  for  calculating  the  Rayleigh  wave  velocity  with  either  shorted  or  free  boundary 
conditions. 

The  shorted  velocity  calculations  assumes  a  massless,  pet  My  conducting  layer  on  the 
surface  of  the  crystal.  The  boundary  condition  which  must  be  satisfied  is  that  <(>=0  at  x} =0.  The  free  surface 
boundary  conditions  are  such  that  the  potential  <f>  and  the  normal  component  of  the  displacement.  D3,  are 
continuous  at  the  surface.  Furthermore,  <t>  must  satisfy  Laplace’s  equation  above  the  surface,  resulting  in 

<t>  =  <t>*  exp(-kxj)  exp  (— ik(Xi  -  Vt)),  (x3  >  0)  <22) 


"‘Surface  Wave  Filters,”  Matthews,  Herbert,  John  Wiley  &  Sons,  New  York  (1977). 


In  both  calculations  the  mechanical  boundary  conditions  are  the  same,  that  there  be  no  force 
component  in  the  x3  direction,  or 

(23) 

Li  =  Tj:  =  Tn  =  0  at  Xi  =  0 

A V/V  for  different  crystal  types  with  various  orientations  were  calculated.  Results  have  been 
found  to  be  in  good  agreement  with  experiments. 

For  example,  on  ST-quartz, 

V.,  =  3.1586  xIO3  m/s 

V,  =  3.1569  x  103  m/s  (24) 

A  V/V  =  5.4  x  10'4 

To  insure  the  suitability  of  the  cuts  described  in  paragraph  (8)  above  for  SAW  applications,  the 
coupling  coefficients  have  been  calculated  for  these  cuts  and  are  summarized  in  Table  3. 

TABLE  3.  PROPAGATION  CHARACTERISTICS  OF  SELECTED  ORIENTATIONS 


& 


Angles  of  ZTCF'11,  degrees 
(S  &  T’s  program) 

Velocity 

(msec) 

K  l 

(x10J) 

Power  FI 

Angle  (Deg 

Phi 

Theta 

Psi 

6 

26 

136.31 

3296.84 

1.12 

-0.3 

6 

27 

135.93 

3293.60 

1.12 

-0.2 

6 

28 

135.59 

3290.63 

1.12 

-0.1 

7 

26 

135.99 

3303.33 

1.12 

-0.5 

7 

27 

135.64 

3299.70 

1.12 

-0.4 

7 

28 

135.27 

3296.33 

1.12 

-0.3 

8 

26 

135.74 

3310.15 

1.12 

-0.7 

8 

27 

135.36 

3306.11 

1.12 

-0.6 

8 

28 

134.97 

3302.32 

1.10 

-0.5 

1 

26 

137.78 

3268.80 

1.10 

+0.7 

1 

27 

137.48 

3267.44 

1.10 

+0.9 

1 

28 

137.17 

3266.36 

1.10 

+1.0 

0 

26 

138.07 

3264.09 

1.12 

+0.9 

0 

27 

137.78 

3263.09 

1.10 

+1.1 

0 

28 

137.49 

3262.35 

1.10 

+1.2 

-1 

26 

138.37 

3259.65 

1.10 

+1.1 

-1 

27 

138.09 

3259.01 

1.10 

-1.3 

-1 

28 

137.80 

3258.64 

1.08 

+1.5 

TABLE  3.  PROPAGATION  CHARACTERISTICS  OF  SELECTED  ORIENTATIONS  (CONT) 


Angles  of  ZTCP",  degrees 
(S  &  T’s  program) 

Velocity 

(msec) 

ie 

(x10‘3) 

Power  Flow 

Angle  (Degrees) 

Phi 

!  Theta 

Psi 

14 

40.195 

3298.60 

0.96 

-7.7 

14 

40.415 

3306.67 

0.96 

-8.1 

14 

41 

40.64 

0.94 

-8.6 

15 

39 

39.79 

wBtfk 

0.96 

-7.8 

15 

40 

40.00 

3310.14 

0.94 

-8.3 

15 

41 

40.23 

3319.09 

0.98 

-8.6 

16 

39 

39.4 

3305.38 

0.96 

-8.0 

16 

40 

39.605 

3314.03 

0.98 

-8.4 

16 

41 

39.825 

3323.15 

0.92 

-9.0 

e.  Power  Flow  Angle 

The  power  flow  angle  for  a  particular  direction  of  propagation  is  an  important  design 
parameter.  While  the  phase  fronts  always  remain  parallel  to  the  source  transducer,  the  wave,  as  a  whole, 
does  not  propagate  perpendicular  to  the  wave  fronts  (see  Figure  7).  This  is  a  characteristic  of  anisotropic 
substrates  where  the  phase  velocity  is  asymmetric  about  the  propagation  direction;  i.e.,  V(i/>  -t-  AiA)  ¥=  \l(ih  - 
A</r).  The  major  problem  which  arises  is  that  the  acoustic  beam  may  steer  off  of  the  desired  propagation 
track,  missing  the  output  trandsucer  unless  it  is  properly  designed. 

The  power  per  unit  width  carried  in  a  surface  wave  is  found  by  integrating  the  mechanical  and 
electrical  Poynting  vectors  to  obtain 


P,  = 


<M)idx3 


,i  =  1.2 


(25) 


Pi  gives  the  power  flow  perpendicular  to  the  wave  front  and  P2  gives  the  power  flow  parallel  to 
the  wave  front.  P5  =  0  for  the  Rayleigh  wave  which  is  confined  to  the  surface.  The  power  flow  angle  may  be 
defined  as 

e  =  arctan(P2/Pi)  »  P2/P.  for  P2 P,  (26) 


The  power  flow  angles  are  calculated  using  either  the  perturbation  programs  or  the  finite 
difference  routine.  Table  4  gives  results  of  calculations  for  the  ST-Cut  quartz.  Note  that  for  <i>  -  0,  the  power 


flow  angle  is  zero,  within  the  single  precision  accuracy  used,  as  a  result  of  the  crystal  symmetry.  Power  flow 
angles  as  high  as  20  degrees  are  not  uncommon  on  quartz. 


TABLE  4.  POWER  FLOW  ANGLES  ST-CUT 


' !> 

Transverse/Incident  Power 

Power  Flow  Angle  e 

0° 

6.159  x  10'5 

0 

10° 

6.145  x  10‘2 

3.5  degrees 

20° 

9.502  x  10‘2 

5.4  degrees 

The  beam  steering  angle  was  calculated  for  selected  cuts  of  quartz  with  TCF  ~  0.  The  analytical 


results  were  compared  with  experimental  measurements.  Table  3  contains  the  results  of  calculations  for  the 
selected  orientations  described  in  paragraph  (8)  above. 


f. 


Proximity  and  Excitation  Strength  of  Bulk  Acoustic  Waves  (BAW)  Spectrum 


The  purpose  of  this  task  is  to  analytically  determine  the  strength  of  spurious  signals  caused  by 
BAW.  With  this  knowledge,  one  can  predict  the  quality  of  the  SAW  device  which  has  been  selected  for  its  zero 
TCD  cut. 


Interdigital  transducers  do  not  couple  all  of  their  field  energy  into  surface  waves.  Bulk  waves 
are  also  generated  at  various  frequencies  in  the  piezoelectric  medium.  These  bulk  waves  can  bounce  off  of 
the  bottom  surface,  undergoing  mode  conversions  in  the  process,  and  be  received  by  the  output  transducer, 
resulting  in  unacceptable  spurious  signals.  In  most  cases,  this  problem  can  be  virtually  eliminated  by  proper 
preparation  of  the  bottom  surface. 

A  transducer  with  periodicity  x  excites  surface  waves  at  center  frequencies 

fsAw  =  (2n  +  1)VSaw/a  (27) 


93 


where  VSA*  is  the  surface  wave  velocity  and  n  is  an  integer,  n = 0, 1 , 2 . The  coupling  to  higher  harmonics 

depends  on  transducer  design.  Bulk  waves  are  radiated  into  the  medium  at  an  angle  &  at  a  center 
fundamental  frequency 

f„  =  V„(0)/[A  cos(0)]  (28) 

where  Vb(0)  is  the  velocity  of  the  bulk  wave  in  that  particular  direction.  Equation  28  is  a  statement  of  the 
Bragg  condition  (see  Figure  8).  If  the  medium  is  isotropic  for  each  mode,  Vb  is  a  constant.  In  the  general 
case,  Vb  depends  on  the  direction  of  propagation  and  hence  on  0.  To  calculate  Vb(0),  the  matrix  of  elastic 
constants  are  rotated  through  the  angle  e  in  the  standard  way.  The  quasi-longitudinal  velocity  and  quasi- 
shear  velocities  are  calculated  as1: 

Vi  =  i[Cn(0)  +  eii2(0)/en(0)]/p}1/2 

Vsh  =  |[C6b(0)  +  eis2(0)/«i  j(0)]/p},/2  (29) 

Vsv  =  f[CAs(0)  +  e  i  s2(0)  / « 1 1(0)]/ p} 1/2 

- - A - ► 


Figure  8.  Bragg  Condition  for  Excitation  of  BAW 
The  inverse  values  of  these  velocities,  plotted  as  a  function  of  0,  form  the  inverse  velocity  curves. 

Inverse  velocity  plots  as  a  function  of  propagation  are  particularly  useful,  for  as  long  as 
Vb/cos  0,  the  effective  velocity  of  the  bulk  wave  on  the  surface,  is  greater  than  Vsaw,  the  effects  of  bulk 
mode  generation  may  be  suppressed  by  appropriate  design  and  by  proper  conditioning  of  the  bottom 
surface.  Should  Vb/ cos  0  be  less  than  VSAw  for  any  bulk  mode,  the  SAW  may  couple  to  the  BAW,  resulting  in  a 
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leaky  surface  wave.  Inverse  velocity  plots  have  been  made  for  orientations  with  promising  SAW 
temperature  characteristics  to  check  for  the  possibility  of  leaky  surface  waves  and  minimum  values  of  fb 
calculated  and  compared  with  f  SAW. 

The  polar  plots  of  the  inverse  velocities  for  a  (YX  wit)  0/27/137.8  and  7 /27/135.59  are  shown  in 
Figure  9  and  Figure  10,  respectively.  The  values  of  (1  /vb»ik)  for  the  two  shear  modes  and  one  longitudinal 
mode  for  different  propagation  directions  (6)  into  the  crystal  are  also  shown  in  Figures  9  and  10.  The 
inverse  surface  wave  velocity  for  0/27/137.8  is  3.06  x  10"\  that  for  7/27/137.8  is  3.03  x  10"\  These  values  are 
larger  than  the  maximum  (1  Nbalk)  cos  e  «2.9  x  IQ"4) .  therefore,  the  analysis  indicates  that  a  leaky  mode  does 
not  exist  (see  Table  3  for  the  SAW  velocities). 

BULK  WAVE  INVERSE  VELOCITIES 

v  in-3  FOR  (YXWLT)  00/  27.00/137.80 


Figure  9.  Polar  Plots  of  Inverse  velocities  for  a  (YXwIt)  0/27/137.8 


BULK  WAVE  INVERSE  VELOCITIES 

x10-3  FOR  (YXWLT)  7.00 /  27.00/135.59 


-0.3  -0.2  -0.1  0.0  0.1  0.2  0.3 

X10'3 

Figure  10.  Polar  Plots  of  Inverse  Velocities  for  a  (YXwlt)  7/27/135.59 


g.  Degeneracies 

Degenerate  waves  occur  when  the  physical  constants  are  such  that  equations  (30)  and  (32) 


decouple. 


r„  -  pV.2 

rI2 
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r.2 
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r„ 

r„ 

r3J  -  pV,: 
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r  54 

r,4* 
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r24 

at 

r«4 

a* 

96 


Tu  —  C55  b  +  2cis  b  +  Cm 

r22  =  C44  b2  +  2c46  b  +  C66 
r»  =  Cjj  b2  +  2cjs  b  +  Css 

r  12  =  C4s  b2  +  (Cm  +  Css)  b  +  Cit 

Tu  =  C35  b2  +  (c  1 3  +  Css)  b  +  Cis 

r2,  =  C34  b2  +  (Css  +  C45)  b+  Cs« 

Ts4  =  — («33  b2  +  2*13  b  +  «ll)  (31) 

Ti4 =  635  b2  +  (61s  +  831)  b  +  On 

r24  =  634  b2  +  (614  +  836)  b  +  61s 

Ts4  =  633  b2  +  (613  +  63s)  b  +  Bn 

The  boundary  conditions  become,  in  matrix  form, 

...(C3311 + Cj3i3 b<m> aim] + (8133+8333 b"*1) «4<b> ...  '  Ci  {32) 

...(Cji,i  +  C„„  b,n”) «,« +(8m  +  833,  b|B>)  as'"' ...  C2 

...(C321I+  C23I3  b'"”)  Ol"*'  +  (8.32+8332^"")  Os*"’ ...  C,  =0 

...(83,,  +e,„  b "*>)  «'"'-(<„  +<3,b|,a*-i<.)«4(m'...  Cs 

-  J  La 

The  definition  of  the  constants  and  derivation  of  the  equations  are  given  in  Appendix  C.  The 
condition  for  decoupling  requires  that  the  matrix  in  equation  30  has  zero  elements  such  that  independent. 
non-Raleigh  wave  solutions  may  exist.  Equation  30  may  decouple  in  many  ways.  If,  for  example,  r,  2  =  Tu  = 
r34  =  0,  u,  and  u3  are  found  to  be  coupled  and  u2  and  <t>  are  found  to  be  coupled,  however  Ui  and  us  are 
decoupled  from  u2  and  <t>.  If  the  physical  constants  are  such  that  these  two  solutions  are  not  coupled  through 
the  boundary  conditions  (equation  32)  then  we  find  the  Rayleigh  like  wave  (Ui  and  u3)  is  not  coupled  to  <t>,  the 
potential  term,  and  cannot  be  excited  by  electrodes  in  this  cut. 

The  displacement  u2  which  is  coupled  to  <t>  is  called  the  Bleustein  —  Gulyaev  wave  and  is 
excited  by  interdigital  transducers.  These  two  waves  (u},  Ui  and  u2,  <t>)  are  degenerate  as  they  propagate  with 
the  same  velocity.  The  simplest  method  to  determine  whether  equation  30  has  decoupled  for  a  particular 
orientation  is  to  calculate  the  matrix  of  equation  30.  This  is  presently  accomplished  using  the  Rayleigh  wave 
velocity  calculation  program,  which  calculates  and  prints  the  matrix  in  equation  30.  Because  of  the  variety 
of  special  cases'  which  may  arise,  some  of  which  may  ultimately  prove  useful,  each  case  in  which  we  find 
the  wave  equation  decouples  will  be  considered  on  an  individual  basis. 

Propagation  Pfezoeier*-'  Surface  Waves  on  Cubic  and  Hexagonal  Crystals,”  Cambell  &  Jones.  J.A.P., 
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The  analytical  results  of  this  section  allow  us  to  determine  which  modes  can  be  generated  on 
the  selected  zero  TCD  cut  of  quartz.  At  the  three  areas  tabulated  in  Table  2,  only  Rayleigh  wave  modes  were 
found. 

h.  Sensitivities  Due  to  Crystal  Misorientation 

In  cutting  quartz  and  aligning  masks  on  it,  there  is  always  some  maximum  achievable 
accuracy.  Thus  it  is  useful  to  know  how  all  of  the  acoustic  quantities  considered  vary  with  angle.  Quantities 
such  as  TCD,  phase  velocity,  power  flow  angle,  aV/V,  bulk  wave  spectrums,  and  bulk  wave  velocity 
surfaces,  are  of  interest  to  this  program.  These  quantities  can  be  accurately  determined  by  directly 
calculating  the  quantities  at  y>  =  (y>„  +  A<#>) ,  ©  =  (0„  +  A0),  and  y- = (y-D +  Ayr)  with  the  same  computer  program 
discussed  in  paragraph  1,  where  <#>„,  0„  and  y-0  are  the  desired  angles,  A y>,  A0  and  A  y-  are  the  actual 
directional  deviation  from  the  desired  direction  due  to  the  fabrication  tolerance.  It  is  important  to  calculate 
the  sensitivity  of  the  parameters  to  the  crystal  misorientation;  i.e.,  amount  of  change  of  a  function  as  a  result 
of  small  angular  misorientation. 

For  case  of  the  power  flow  angle  (PFA) 

d[PFA  (y>  =  0,  0  =  42.75,  *  =  10°)]/dyr  f33* 

s  (PFA  (yr  =  20°)  -  PFA  (*  =  0°))/20° 

s  0.27°/degree. 

The  quantities d[PFA(y.=0°,  0=42.75,  *>=10°)]/dy> and d[PFA(y>=0°,  0  =  42.75°,  ^=1O°)]/d0are computed  in 
an  equivalent  manner. 

In  the  case  of  TCD: 

d(TCD)  =  +  JOCQId*  +  8QCBLd#  (34) 

ou  o<t>  oyi 

=  [TCD(y>„,  0O,  yg  -TCD(*0, 0„+ A0,  *.)]/ A« 

96  (35) 

=  (TCD(0o,  0o,  yro)-TCD(<^o+ A0, 0„,  y».)]/Ay» 

o<J> 

=  [TCD(y>o,  00.  y»0)  -  TCD(y>0, 0„,  y»0 + Ayr)]/  Ay- 

All  the  TCD  (y>,  0,  y -)  and  TCD  (y>,  0,  y-  +  Ay-)  etc.,  are  calculated  with  the  computer  programs  discussed  in 
paragraph  1.  Therefore,  all  the  9(TCD)/30, 3(TCD)/ayr,  3(TCD)/30  can  be  accurately  calculated.  Hence  the 
sensitivities  due  to  crystal  misorientation  are  determined. 

Calculation  of  the  angular  dependence  on  the  first,  second,  and  third  order  TCDs  is.  of  course, 
our  primary  task.  Of  these  three  quantities,  the  first  order  TCF  is  most  sensitive  to  angular  variation.  The 


angular  dependence  on  these  parameters  were  calculated  in  the  same  way  as  all  of  the  other  quantities  but 
on  a  much  smaller  angular  grid  (Ad,  a<*>  and  a< p)  about  the  zero  TCD  locus.  The  size  of  this  grid  varies, 
depending  on  the  magnitude  and  smoothness  of  the  variation  about  each  point  on  the  locus  of  zero  first 
order  TCD  orientations  as  verified  by  experiment  and  calculation.  In  case  the  function  is  smooth  or  linearly 
varying  with  angles,  large  a  angles  may  be  used. 

By  performing  the  above  mentioned  calculations,  we  have  precise  information  on  the 
sensitivities  due  to  crystal  misorientation.  This  information  will  allow  us  to  impose  a  practical  tolerance 
limit  on  fabrication  and  still  be  able  to  achieve  the  required  superior  performance  specification. 

Quantities  such  as  velocity  (Table  3),  power  flow  angles  (Table  3),  BAW  spectrum  (Figures  50 
and  51),  coupling  coefficients  (T able  3),  and  second  and  third  order  TCFs  (Table  2)  do  not  vary  quickly  with 
angle.  This  is  not  the  case  for  TCF1".  Table  5  contains  a  summary  of  aTCF"’/a</f.  The  large  values  of 
3TCF"'/a*  impose  strict  fabrication  tolerances  on  the  SAW  cuts  and  mask  alignment.  Fabrication  accuracy 
to  within  6  minutes  is  required  to  keep  the  total  temperature  variation  due  to  TCF111  within  45  ppm  for 
aTCF"7a«A  =  3(PPM/°C)/degree  over  the  temperature  range  -50°C  to  100°C.  Table  6  contains  summaries  of 
3TCF"730  and  aTCF"7 do.  These  values  impose  fabrication  tolerances  on  the  rotated  quartz  plate  angles  <t> 
and  e  of  12  minutes  to  keep  the  total  temperature  variation  due  to  aTCF'"(15/40/40)/a<<>  within  45  ppm  over 
the  temperature  range  -50°C  to  100°C.  This  linear  temperature  variation  may  be  compensated  for  by  varying 
ib  on  any  particular  cut  if  all  other  cut  parameters  vary  slowly  with  angle. 


TABLE  5.  aTCF“»/a^  FOR  SELECTED  CUTS 


Angles  of  ZTCF"1,  degrees 
(S  &  T's  program) 

aTCF"7a* 

Phi 

Theta 

Psi 

6 

26 

136.31 

+2.7  (ppm/C°)/degree 

6 

27 

135.93 

+2.7 

6 

28 

135.59 

+2.7 

7 

26 

135.99 

+2.7 

7 

27 

135.64 

+2.7 

7 

28 

135.27 

+2.7 

8 

26 

135.74 

+2.7 

8 

27 

135.36 

+2.7 

8 

28 

134.97 

+2.7 

1 

26 

137.78 

+2.8 

TABLE  5.  dTCF"7a*  FOR  SELEC1 

[ED  CUTS  (C0NT) 

Angles  of  ZTCF"1,  degrees 
(S  &  T's  program) 

aTCF'"/0(/, 

Phi 

Theta 

Psi 

1 

27 

137.48 

+2.8 

1 

28 

137.17 

+2.8 

0 

26 

138.07 

+3.0 

0 

27 

137.78 

+3.0 

0 

28 

137.49 

+3.0 

-1 

26 

138.37 

+3.0 

-1 

27 

138.09 

+3.0 

-1 

28 

137.80 

+3.0 

14 

39 

40.195 

-3.5 

14 

40 

40.415 

-3.5 

14 

41 

40.64 

-3.5 

15 

39 

39.79 

-3.5 

15 

40 

40 

-3.5 

15 

41 

40.23 

-3.5 

16 

39 

39.4 

-3.7 

16 

40 

39.605 

-3.7 

16 

41 

39.825 

-3.7 

TABLE  6.  aTCF"7a*  AND  aTCF"7a0  FOR  SELECTED  CUTS 


Angles  of  ZTCF'” 

(S&T's  Program),  Degrees 

3TCF"'/a<i> 

aTCF"7ao 

Phi 

Theta 

Psi 

7 

m 

135.64 

-0.7(ppm/C°)/degree 

-0.5(ppm/C°) /degree 

0 

137.78 

-0.8 

-0.8 

15 

40 

40.00 

+1.5 

-0.7 

2.  WAFER  FABRICATION 


During  the  developmental  phase  of  the  doubly  rotated  cut  SAW  Devices  program,  a  large  number  of 
wafers  with  30  different  crystal  orientations  were  fabricated  and  evaluated.  To  perform  this  task  with  good 
quality  control,  minimum  cost  and  in  a  short  time,  Motorola  developed  the  techniques  and  processes  to 
fabricate  the  doubly  rotated  cut  of  quartz  wafers  internally.  The  crystal  boules  are  supplied  by  Motorola. 
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Carlisle,  Pennsylvania;  the  X-ray  wafer  cutting  is  performed  at  Motorola’s  Semiconductor  Group;  polishing, 
fabrication  and  testing  are  performed  at  the  Motorola  Government  Electronics  Division.  The  methods  used  to 
fabricate  the  quartz  wafers  are  described  in  the  following  sections. 

A  computer  program  was  developed  to  calculate  the  incident  angles  and  reflected  angles  for  any 
given  cut  of  crystal.  The  basic  mathematical  relations  were  derived  in  R.A.  Heising's  “Quartz  Crystals  for 
Electrical  Circuits."  The  relations  were  formulated  to  computer  programs.  The  flow  diagram  is  shown  in 
Figure  1 1 .  The  results  are  shown  in  the  printout  (see  Appendix  F).  The  incident  angle  G,  exit  angle  (G + G’),  are 
defined  for  each  plate  position.  The  plate  positions  are  determined  by  the  direction  of  the  rotated  axis,  (±1, 
±2,  ±3),  relative  to  the  measuring  stage.  Once  the  angles  of  (YX  wit)  <t>lol<ii  are  defined,  all  of  the  reflection 
angles  will  be  calculated  for  the  different  reflection  planes,  defined  by  Miller  indicies.  The  useful  reflection 
planes  that  provide  intense  reflection  with  low  skewed  angles  (<5  degrees)  are  selected  for  printout.  These 
results,  after  modification  due  to  change  in  standards,  are  used  to  check  the  accuracy  of  the  cuts. 

a.  Face  Definitions 

The  following  is  an  outline  of  the  procedure  to  be  used  in  this  program  to  make  a  doubly 
rotated  cut  (YX  wit)  4>lel 0.  The  quartz  bars  have  four  lumbered  faces  with  the  minus  X  axis  marked.  The 
opposite  face  is  marked  by  coloring  it  with  a  magic  marker  (blue).  One  other  face  not  opposite  the  -X  face  is 
colored  also  with  a  magic  marker  (red).  The  red  face  is  now  defined  to  be  the  +Z  axis  as  in  Figure  12.  The 
direction  of  the  +Y  axis  may  be  found  by  using  the  right-hand  rule  (+X  crossed  into  +Y  gives  +Z).  The  +Y  axis 
will  lie  along  the  length  of  the  crystal.  Note  that  there  are  two  ways  to  set  up  the  axes  on  the  crystal 
corresponding  to  choice  of  the  rad  face  (see  Figure  12). 

b.  Running  the  X-Ray  Program 

The  X-ray  orientation  program  XRAY  is  run  with  the  angles  <*>/0/O.  The  sequence  of 
instructions  on  the  Honeywell  560  under  the  CP-V  operating  system  is  as  follows: 

'.SET  F:103  DC/MILDAT  .538;  IN 

(This  instruction  assigns  to  unit  103  the  tile  containing  all  of  the  Miller  indicies  to  be  searched.) 

ISET  F:104  DC/ERRDAT  ,538;IN 

(This  instruction  assigns  to  unit  104  the  file  containing  all  of  the  angle  perturbations  defined  in 
Heising  for  use  with  the  Laue  photographs.) 

IXRAYL  M.538 

(This  instruction  loads  and  begins  execution  of  the  program.) 

The  program  will  ask  for  the  angles  <t>lolO.  These  angles  are  entered  in  3G  format  [i.e.,  10.0, 3.0, 
0.0  (carriage  return)]  and  the  results  of  the  X-ray  analysis  printed. 
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CALCULATE  ROTATION 
MATRIX  FOR  (YXWLT)4>/0/^ 


Figure  11.  Flow  Diagram  of  Calculation  for  Reflected  Angles  in  Doubly  Rotated  Cut  of  Quartz 


I 
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2510-5 


Figure  12.  Two  Possible  Ways  of  Defining  the  Crystal  Axes  on  a  Lumbered  Bar 


c.  Analysis  of  X-Ray  Program 

The  X-ray  program  used  follows  the  development  of  Heising’s  which  uses  a  book  written 
before  the  1949  IRE  Standard  we  use  now.  His  X.  Y,  and  Z  axes  will  be  denoted  as  X,  Y,  and  Z.  The  relation  to 
our  standard  is  as  follows: 


Heising 

1949  IRE  Standard 

+X 

-X 

+Y 

-Y 

+Z 

+Z 

Furthermore,  Heising’s  incident  beam  comes  from  the  left,  while  our  incident  beam  comes  from  the  right. 
When  these  differences  are  accounted  for,  the  "position  chart"  shown  in  Figure  13  must  be  used  instead  of 
the  charts  Heising  uses  in  order  that  our  X-ray  machine  and  the  1949  standard  may  be  used.  P,,  P2,  and  P3  in 
the  "position  chart”  are  the  +X,  +Y,  and  +Z  axes  of  the  rotated  plate.  An  important  note  of  caution:  the  <t> 
crystal  face  off  of  which  we  reflect  the  X-ray  beam  has  Y  or  -Y  as  its  normal.  Since  we  may  not  assume  that 
the  crystal  faces  are  parallel,  we  must  propagate  the  surface  wave  on  this  face.  This  is  important,  because  to 
achieve  the  same  crystal  properties,  we  must  rotate  by  about  the  +Y  axis,  resulting  in  a  change  of  the  sign 
of  the  third  rotation  (see  Figure  14). 
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INCIDENT  BEAM 


P3 


POSITION  1 


POSITION  3 


POSITION  4 


Pi,  P2,  and  P,  are  the  rotated  plate  axes  X,  Y,  and  Z. 


2510-7 


Figure  13.  X-ray  Position  Chart 


THIRD  ANGLE  ROTATION 


iv 


If  the  angles  a  and  p  in  Figure  16  are  measured  using  the  (003)  plane,  which  is  parallel  to  the  X  and  Y  axes, 
is  given  by: 


A\l/  =  -tan  '[(-cos  <t>  sin  a  -  sin  <t>  sin  /3)/(sin  e  sin  0  sin  a  -  sin  0  cos  <t>  sin  p  +  cos  0  >/l  -sin2«  -  sin20)] 

(36) 


Thus,  given  a  ^  desired  which  we  wish  to  obtain,  we  must  actually  rotate  the  mask  by  =  ii/6a -a*  as 
shown  in  Figure  17. 

e.  Determination  of  a  and  p 

The  most  appropriate  crystal  plane  to  use  for  determining  a  and  p  is  the  003  plane.  This  plane 
has  all  of  its  atoms  lying  in  planes  perpendicular  to  the  Z  axis,  and  for  a  perfect  Z  face  the  X-ray  deflection 
angles  will  be: 

G  =  Angle  of  incidence  =  25  degrees  19  minutes 
G  +  G’  =  Exit  angle  =  50  degrees  38  minutes 

For  all  positions  on  the  X-ray  machine  with  the  +Z  or  -Z  face  being  X-rayed.  To  determine  a  and  /3for  the  -Z 
face,  we  use  the  following  relationships. 

a  =  1/2  (G  -X  TOWARD  INCIDENT  “  G*X  TOWARD  INCIDENT)  (37) 

P  ~  1',.  (G*y  toward  incident  “  G-y  toward  incident)  (38) 


TRUE  INTERSECTION 


To  determine  a  and  p  for  the  +Z  face  we  use 

at  =  M2  (G.X  TO  INCIDENT  —  G-X  TO  INCIDENT)  (39) 

0  -  M2  (G  -Y  TO  INCIDENT  ~  G+Y  TO  INCIDENT)  (40) 

a  and  0  are  the  a  and  p  we  would  have  measured  for  a  -Z  face  that  is  parallel  to  the  +Z  face  measured, 
which  means  that  equation  (36)  is  still  valid  for  a  and  0. 

The  quartz  crystal  will  now  have  the  orientation  (YX  wit)  <plol  -  A</»(a,  p)  if  we  use  the  -Z  edge 

as  a  reference  and  (YX  wit)  <t>lol  -  A0(a\  0)  if  we  use  the  +Z  edge  as  a  reference.  The  computer 

X-ray  program  for  these  angular  orientations  must  be  used  if  we  are  to  get  an  accurate  measurement.  The 
integrity  of  these  measurements  of  a  and  p  is  checked  by  comparing  1 12  (G-x  +  G*x)  and  1 12  (G-v + G.v)  with 
the  theoretical  result,  GAvc  =  25  degrees  19  minutes. 

f.  Procedure  Used  to  Zero  X-Ray  Machine: 

(1)  Use  reference  quartz  plate  marked  10.1 
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(2)  Set  X-ray  exit  angle  to  26.5  degrees  (marked  AT  cut  on  machine) 

(3)  Set  plate  in  X-ray  machine  with  marking  toward  you  (X-ray  bounces  off  unmarked 
surface)  and  blue  line  up  (width  direction  vertical).  (See  Figure  18). 

(4)  Adjust  with  clutch  till  X-ray  reading  occurs  at  incident  angle  13  degrees  18  minutes 

(5)  Set  plate  in  X-ray  machine  with  marking  toward  you  (X-ray  bounces  off  unmarked 
surface)  and  blue  line  down  (width  direction  vertical).  (See  Figure  19). 

(6)  Angle  of  incidence  should  '•••  •.*  at  13  degrees  18  minutes  +  a  angle. 

(7)  Adjust  clutch  until  jspgM  ^id'nce  is  at  13  degrees  18  minutes  +  1/2  A  angle. 

(8)  Check  alignment  by  g  in  positions  of  steps  3  and  5.  The  average  of  the  two 

readings  should  be  13  degrees  18  minute'  A:<  s  ■*  screpancy  represents  error  in  zeroing  the  machine  and  may 
be  made  less  than  iess  than  1  minute. 

(9)  If  the  10.1  cut  is  not  available,  the  same  procedure  may  be  used  if,  for  the  cut  to  be  used, 
the  theoretical  angles  of  incidence  are  equal  for  the  particular  crystal  plane  used  in  both  positions. 


BLUE  LINE 


Figure  18.  Position  of  Step  3, 10.1  Alignment  Wafer 


Figure  19.  Position  of  Step  5, 10.1  Alignment  Wafer  Rotated  by  180  Degrees  About  -Y  Axis 


******&* 
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g.  Cutting  Techniques  for  Doubly  Rotated  Orientations 


A  slicing  machine  built  by  Meyer  and  Burger  Company,  model  TS3,  was  used  to  slice  the  quartz 
boule.  The  6  inch  diamond  impregnated  blade  is  supplied  by  Maurice  Dessau,  New  York.  The  drift  accuracy 
of  the  2  inch  cut  is  approximately  3  mils  (~5  minutes). 

The  stage  of  the  saw  has  the  capability  to  rotate  in  two  dimensions;  it  is  designed  to  make  the 
doubly  rotated  cut.  The  adjustment  accuracies  are  graded  to  0.5  minute.  The  X-ray  beam  resolution  is  about 
2  minutes. 

The  combined  accuracy  of  the  wafers  obtained  experimentally  has  been  better  than  15 
minutes.  Iterative  adjustment  of  the  cutting  stage  can  bring  the  cutting  accuracy  to  within  10  minutes.  The 
cutting  procedures  are  the  following. 

(1)  Heat  brass  plate,  glass  plate  and  quartz  boule. 

(2)  Melt  wax  on  brass  plate  and  mount  glass  plate. 

(3)  Melt  wax  on  glass  and  mount  quartz  boule  on  glass. 

(4)  Let  cool  to  room  temperature. 

(5)  Screw  mount  brass  plate  to  slice  machine. 

(6)  Initial  cut  of  boule  along  the  XZ  plane  for  directly  to  the  desired  angle  if  the  boule  was 

lumbered). 

(7)  Determine  true  atomic  plane  using  diffractometer  and  X-ray  computer  program. 

(8)  Adjust  stage  to  correct  for  true  atomic  planes,  and  obtain  the  desired  rotation. 

(9)  Perform  doubly  rotated  cut  of  quartz. 

(10)  Measure  the  reflection  angle  and  check  against  the  calculated  result  of  the  X-ray 
program.  Iterative  adjustment  of  stage  can  be  done  if  cut  accuracy  is  critical. 

These  procedures  are  presently  used  to  cut  the  wafers  with  low  TCF  orientation  required  for 

this  program. 

h.  Wafer  Polishing 

The  polishing  process  for  the  quartz  wafers  was  established  during  this  period.  The  procedure 
includes  the  following; 

(1)  The  crystal  axis  orientation  is  marked  on  the  back  surface  of  the  wafers  with  marking  ink 
and  baked  dry.  The  markings  have  been  demonstrated  to  stay  on  the  crystal  throughout  the 
photolithographic  process  for  easy  identification. 


(2)  The  corners  of  the  wafers  are  ground  prior  to  lapping  and  polishing  to  prevent  corner 


breakage. 


(3)  The  polishing  is  a  two  step  process.  The  first  step  is  a  15  urn  lapping  to  obtain  surface 
flatness  with  Microgrit,  type  WCA,  Size  15,  on  a  LAPMASTER,  Model  24.  The  second  step  is  to  polish  the 
surface  with  Cerium  oxide  on  a  LogiTech,  Ltd.  Model  PM2  polisher.  The  polish  wheel  is  made  of  pitch  wax 
supplied  by  Hacker  Instruments,  Inc.  The  polished  wafers  have  a  surface  finish  with  no  observable  grains  or 
pits  under  a  50X  microscope.  The  wafers  are  suitable  for  surface  wave  application.  X-ray  orientation  was 
performed  before  and  after  the  lapping  and  polishing  processes.  The  change  is  normally  less  than  6  minutes. 
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SECTION  III 

TECHNICAL  DISCUSSION  OF  TASK  II 

1.  INTRODUCTION 

The  search  for  a  temperature  stable  cut  of  quartz  for  application  to  SAW  devices  has  led  to  the 
investigation  of  the  doubly  rotated  cuts.  Theoretical  studies  have  indicated  that  doubly  rotated  cuts  of 
quartz  promise  much  better  temperature  stability  than  the  commonly  used  ST  cut.  Task !  of  this  program, 
which  encompassed  the  first  iteration  calculations  of  the  doubly  rotated  cuts  of  quartz  was  successfully 
completed.  Task  II,  with  major  emphasis  on  experimental  work,  is  discussed  in  this  section. 

In  Task  I  theoretical  analyses  have  been  performed  and  angular  rotations  promising  very  low  TCF'n 
and  TCP2’  have  been  plotted.  Important  SAW  device  design  parameters,  such  as  coupling  coefficient, 
velocity  and  power  flow  angle,  have  also  been  computed  to  characterize  each  area.  As  part  of  Task  II, 
experimental  results  establishing  the  degree  of  correlation  with  theory  have  been  obtained. 

Theoretical  propagation  characteristics,  as  discussed  in  paragraph  2  below,  impose  strict 
fabrication  tolerances  on  the  SAW  cuts  and  mask  alignment  due  to  large  values  of  STCF1  /9i/».  Thus  a  mask 
had  to  be  designed  to  compensate  for  fabrication  errors.  This  design  work  is  presented  in  paragraph  3.  Since 
experimental  technique  is  an  important  criteria  in  the  determination  of  correlation  between  theory  and 
experiment,  all  procedures  followed  as  well  as  equipment  used  are  discussed  in  paragraph  4.  Paragraph  5 
presents  experimental  determination  of  propagation  characteristics  which  illustrate  an  excellent  agreement 
between  theoretical  calculation  and  experimental  results.  Paragraph  6  contains  the  results  of  the  doubly 
rotated  cut  TCF  measurements  made  to  date. 

Theoretical  calculations  have  been  in  good  agreement  with  experimental  results.  Doubly  rotated  cuts 
of  quartz  with  an  improvement  of  TCFt2'  by  at  least  a  factor  of  two  over  the  ST  cut  have  been  obtained.  A 
further  improvement  was  expected  after  a  second  iteration  and  these  results  are  discussed  in  Section  IV. 

2.  THEORETICAL  PROPAGATION  CHARACTERISTICS 

In  cutting  quartz  and  aligning  masks  on  it,  there  is  always  some  maximum  achievable  accuracy.  It  is 
also  useful  to  know  how  all  of  the  acoustic  quantities  considered  vary  with  angle.  Quantities  such  as  TCF, 
phase  velocity,  power  flow  angle,  aV/V,  and  bulk  wave  velocity  surfaces  are  of  interest  to  this  program. 
These  quantities  can  be  accurately  determined  by  directly  calculating  the  quantities  at  <*>  =  (</>0+A </>),  e  = 
(00+Afl),  and  < h  =  (i/r0+A<|r).  Calculation  of  the  angular  dependence  on  the  first,  second,  and  third  order  TCF's 
is,  of  course,  our  primary  task.  Of  these  three  quantities,  the  first  order  TCF  is  most  sensitive  to  angular 
variation  (refer  to  Tables  2, 3  and  5).  Quantities  such  as  velocity  (Table  3),  power  flow  angles  (Table  3), 
coupling  coefficients  (Table  5).  and  second  and  third  order  TCF’s  (Table  2)  do  not  vary  quickly  with  angle. 
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This  is  not  the  case  for  TCF.1'1  Table  5  contains  a  summary  of  3TCF‘/3(/<.  The  large  values  of  3TCF"73i/< 
impose  strict  fabrication  tolerances  on  the  SAW  cuts  and  mask  alignment.  Therefore  it  is  essential  to  design 
a  mask  with  reference  registration  marks  to  accurately  determine  the  transducer  orientation  relative  to  the 
crystal  edges.  These  reference  markers  are  fabricated  on  all  of  the  measured  devices,  so  that  propagation 
direction  is  accurately  determined  to  within  ±25  minutes.  Fabrication  accuracy  to  within  6  minutes  is 
required  to  keep  the  total  temperature  variation  due  to  TCP"  within  45  ppm  for  3TCP"/3^  =  0.3 
(PPM/°C)/degree  over  the  temperature  range  of  -50°C  to  100°C.  Table  10  contains  summaries  of  3TCF"/3<*> 
and  3TCF'"/30.  These  values  impose  fabrication  tolerances  on  the  rotated  quartz  plate  angles  <#> and  e  of  Ac*> 
and  A0  less  than  12  minutes  to  keep  the  total  temperature  variation  due  to  3TCF"  (15/40/40)/3<*>  within  45 
ppm  over  the  temperature  range  of  -50°C  to  100°C  for  example.  This  linear  temperature  variation  may  be 
compensated  for  by  varying  <i>  on  any  particular  cut  if  all  other  cut  parameters  vary  slowly  with  angle.  To 
date  all  other  cut  parameters  have  been  found  to  vary  slowly  with  Phi,  Theta  and  Psi. 

3.  MASK  DESIGN 

A  mask  was  designed  to  take  into  account  the  sensitivity  of  TCP"  due  to  small  variations  in  the  cut 
angles.  The  design  incorporates  rotated  structures.  Each  device  is  offset  with  respect  to  its  neighbor  by  0.2°. 
Three  individual  devices  are  illustrated  in  Figure  20.  The  device  specifications  are  as  follows: 

a.  Transducer  periodicity:  12.192  /xm  (center  frequency  ~260  MHz;  varies  with  crystal 
orientation) 

b.  Delay  time:  360x  (~1.4  /xs)  (varies  with  crystal  orientation) 

c.  Number  of  sets  in  output  transducer:  15 

d.  Electrode  pairs  per  set:  2.25 

a.  Electrode  pairs  in  input  transducer:  24 

f.  Aperture  width:  Input:  70  wavelengths 
Output:  50  wavelengths 


Figure  21  is  a  photograph  ot  a  wafer  fabricated  using  the  above-mentioned  design.  Twelve  oscillators 
propagating  in  different  i fr  directions  are  fabricated  on  a  single  wafer.  The  principal  advantage  of  this  design 
is  the  ability  to  compensate  for  fabrication  errors .  Fi  om  the  experimental  results  plotted  in  Figure  22,  we  are 
able  to  observe  significant  frequency  shifts  due  to  small  variations  in  PSI  angle.  These  results  confirm  our 
theoretical  calculations. 

4.  EXPERIMENTAL  MEASUREMENT  OF  TCF 

Wafers  with  orientations  that  provide  low  temperature  coefficients  of  frequency  were  used  to 
fabricate  SAW  oscillator  devices.  Considerable  care  has  been  taken  to  minimize  fabrication  tolerances.  For 
angles  <t>  and  e,  the  estimated  accuracy  is  within  ±4  minutes;  for  angle  < i>,  the  accuracy  is  within  ±25  minutes. 

The  delay  line  oscillators  described  previously  were  used  to  measure  the  frequency  stability  at 
different  temperatures.  The  experimental  apparatus  is  shown  in  Figure  23.  No  coils  were  used  to  match  the 
devices  in  order  to  eliminate  inductance  changes  in  the  matching  circuit  over  the  temperature  range  tested. 

The  switches  are  designed  to  test  ten  oscillators  in  the  same  temperature  chamber.  The  phase 
stability  of  the  switches  is  essential  to  the  measurements.  The  circuits  are  shown  in  Figure  24.  The  phase 
stability  of  the  switches  was  found  to  be  adequate  for  measuring  the  temperature  stability  of  the  device. 

A  digital  thermocouple  (Fluke  2160A-T)  was  taped  :u  the  bottom  of  the  fixture  to  measure  device 
temperature.  The  Fluke  2160A-T  is  accurate  to  within  ±2°C  over  the  temperature  range  of  -75°C  to  +150°C.  A 
thermometer  was  also  used  to  measure  the  oven  air  temperature.  Semi-rigid  cable  constitute  ?!! 
interconnections.  This  reduces  loss,  shields  against  feedthrough,  and  makes  the  apparatus  less  sot$*v*  s  to 
the  testing  environment.  ANZAC  DS109  power  splitters  and  AVETEK  AWL500M  amplifiers  were  used  in  the 
feedback  loop.  A  Systron  Donner  (PLS  50-1 )  provided  the  DC  power  supply  voltage.  The  mean  supply  voltage 
was  maintained  at  15.000  ±  0.001  volts  during  the  measurements. 

Frequency  measurements  were  taken  every  10°C,  spanning  the  range  from  -55°C  to  +135°C. 
Stabilization  of  temperature  and  frequency  was  attained  for  each  measurement  before  data  was  taken.  This 
ensures  that  the  device  is  in  thermal  equilibrium  with  its  environment.  The  total  experimental  error  is 
estimated  to  be  less  than  ±10  PPM. 

5.  EXPERIMENTAL  PROPAGATION  CHARACTERISTICS 

The  propagation  characteristics  of  selected  devices  were  measured  to  verify  the  calculations.  The 
experimental  measured  velocity  is 

v  =  f  k  A  =  12.192  fim 

Table  7  lists  the  experimental  velocities  for  the  various  cuts.  The  deviation  (up  to  0.2%)  is  in  part  caused  by 
the  slight  deviation  of  crystal  orientation  due  to  the  fabrication  process,  and  in  part  caused  by  the 
uncertainty  of  center  frequency  due  to  the  unknown  phase  shift  in  the  feedback  loop  (±0.25  MHz). 
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Figure  21.  Photograph  of  a  Fabricated  Wafer  With  Propagation  Directions  ol 
8°277270547133°54'  ±  n(0.2°) 


Frequency  Temperature  Curve  With  Difterent  Angles 


POWER 


Figure  23.  Measurement  System 
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The  power  flow  angle  of  the  doubly  rotated  cut  quartz  wafer  has  been  measured  by  the  laser  probe 
technique.  The  measured  result  is  in  good  agreement  with  calculations. 

TABLE  7.  EXPERIMENTAL  VELOCITIES  OF  CUTS 


Phi 

Theta 

Psi 

Measured 

Velocity 

Calculated 

Velocity 

-1.05 

28.0667 

136.534 

3257.4 

3260.9 

-0.9667 

26.233 

138.449 

3256.8 

.  3259.4 

-0.133 

28.1 

137.692 

3259.0 

3261.5 

-0.033 

26.7 

138.859 

3262.1 

3262.5 

0.633 

26.15 

137.016 

3267.6 

3267.7 

5.583 

27.833 

135.194 

3289.7 

3288.9 

5.583 

27.833 

135.994 

3290.7 

3288.9 

6.0 

26.9667 

135.812 

3288.9 

3293.7 

6.067 

25.933 

133.099 

3298.8 

3299.4 

7.41 

27.83 

134.2 

3296.8 

3299.1 

8.033 

26.9667 

134.618 

3304.1 

3306.1 

14.2833 

39.1167 

40.227 

3294.8 

3301.4 

14.2833 

39.1167 

40.627 

3296.8 

3304.6 

15.25 

39.2 

39.6187 

3300.7 

3303.4 

15.3 

40.6833 

40.0308 

3314.0 

3317.3 

The  block  diagram  of  the  experimental  set-up  is  shown  in  Figure  25. 


In  this  procedure,  the  first  order  deflecting  light  due  to  the  presence  of  the  acoustic  wave  is  measured  by  the 
photomultiplier.  The  light  intensity  is  proportional  to  the  acoustic  power,  while  the  angle  of  deflection  is 
given  by: 

Sin  ft,  =  Sinft,  +  ^r 

Where  ft  =  angle  of  deflection  of  nth  order 
ft  =  angle  of  specular  reflection 
x  =  optical  wavelength 
a  =  acoustic  wavelength 
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The  angular  relation  is  demonstrated  in  Figure  26.  The  optical  beam  is  provided  by  a  He-Ne  laser  with  spot 
size  of  -100/u.  The  acoustic  waves  are  generated  by  a  260  MHz  transducer  with  a  equal  to  12.192m.  During 
the  measurement,  the  SAW  device  is  translated  up  and  down  so  that  the  optical  beam  is  scanned  across  the 
acoustic  path  to  detect  the  acoustic  beam  intensity  distribution.  The  test  equipment  is  shown  in  Figure  27. 

Figures  28  and  29  show  examples  of  plots  of  the  relative  acoustic  beam  intensities  indicated  by  the 
deflected  laser  beam.  These  plots  were  taken  in  the  near-  and  far-field  regions,  respectively.  The  exact 
distance  between  the  near-field  scan  and  the  far-field  scan  is  8.131  mm.  The  center  of  the  beams  is  estimated 
to  be  0.05  mm  in  separation. 

The  tTieasured  power  flow  angle  for  this  cut  is 

0P  =  tan'1  -  =  -0.352° 

The  calculated  power  flow  angle  for  this  cut  is 

o'P  =  -0.3° 

The  good  agreement  between  the  measurement  and  the  calculated  results  indicates  that  the  calculation  is 
accurate.  Acoustic  dispersion  and  loss  can  also  be  calculated  from  this  data. 

6  EXPERIMENTAL  MEASUREMENT  OF  FREQUENCY  VERSUS  TEMPERATURE 

A  few  representative  frequency-temperature  measurements  are  presented  in  Figures  30  though  34. 
The  stars  represent  experimental  data  points.  The  solid  lines  are  linearly  regressed  curves  used  to  define  the 
measured  first,  second,  and  third  order  TCF's  for  these  cuts,  given  in  Table  8.  Figures  30  through  34  are 
representative  of  cuts  in  region  (YX  wit)  7/27/135.  Cut  (YX/wIt)  6.57/26.88/134.9  of  Figure  30  has  a  small 
linear  frequency  term  at  room  temperature  and  is  well  suited  for  use  at  both  high  and  low  temperatures.  Its 
second  order  TCF  is  in  good  agreement  with  the  computer  calculations  and  is  considerably  smaller  than  that 
of  ST-cut  quartz  (see  Table  8).  Cut  (YX  wit)  5.58/27.83/135.1  of  Figure  31  displays  a  larger  total  frequency 
variation  over  the  range  shown  but  is  much  more  stable  at  higher  temperatures.  This  illustrates  how  a 
simple  change  of  crystal  orientation  can  be  used  to  temperature-compensate  doubly  rotated  cut  SAW 
devices  for  different  mean  operating  temperatures.  A  slight  rotation  of  1 1>,  as  shown  in  Figure  22,  could  be 
used  to  set  the  first  order  TCF  to  zero  while  slightly  altering  the  second  order  TCF  of  the  device.  This  cut  also 
represents  a  substantial  improvement  over  the  ST-cut  (see  Table  8).  Figures  30  through  34  summarize  some 
of  our  typical  measurements  to  date.  Results  in  this  area  are  in  excellent  agreement  with  the  theory.  One  of 
these  cuts  (8.033/26.967/134.6)  has  temperature  stability  of  approximately  40  ppm  from  0°C  to  130°C;  it  is 
suitable  for  systems  or  weapons  operating  in  elevated  temperatures  (see  Figure  33). 

On  all  cuts  tested,  the  agreement  between  the  experimental  and  calculated  results  has  been  excellent. 
These  results  establish  a  firm  basis  for  performing  the  second  iteration  search  for  the  optimum  TCF 
orientations. 


SIGNAL 

GENERATOR 


AMPLIFIERS 


ATTENUATOR 


SAW 

DEVICE 


FREQUENCY 

COUNTER 


Figure  25.  Experimental  Set-Up  For  Laser  Probe 


Figure  26.  Angular  Relation  of  Incidence  and  Reflected  Beam 


!7.  Test  Set-up  for  Laser  Probing  of  Acoustic  Wave  Beam-Steering  Characteristics 
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RELATIVE  ACOUSTIC  BEAM  INTENSITY 
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Figure  28.  Intensity  Distribution  at  Near  Field 
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RELATIVE  ACOUSTIC  BEAM  INTENSITY 
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Figure  29.  Intensity  Distribution  at  Far  Field 
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EXP  ANGLES  ARE  6.567,  26.883, 134.9 


Figure  30.  Experimental  Frequency  Response  (Test  No,  1) 
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Figure  31.  Experimental  Frequency  Response  (Test  No.  2) 
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Figure  33.  Experimental  Frequency  Response  (Test  No.  4) 
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SECTION  IV 

TECHNICAL  DISCUSSION  OF  2ND  ITERATION 


1.  INTRODUCTION 

The  search  for  a  temperature  stable  cut  of  quartz  for  application  to  SAW  devices  has  led  to  the 
investigation  of  the  doubly  rotated  cuts.  Theoretical  studies  have  indicated  that  doubly  rotated  cuts  of 
quartz  promise  much  better  temperature  stability  than  the  commonly  used  ST  cut.  The  first  iteration  of  tins 
program,  which  encompassed  the  Task  I  calculations  of  the  doubly  rotated  cuts  of  quartz,  and  the  Task  II 
experimental  work,  were  discussed  in  Sections  II  and  III.  The  second  iteration  calculations  and  analysis  and 
experimental  results  are  shown  in  this  section. 

In  Task  I,  theoretical  analyses  were  performed  and  angular  rotations  promising  very  low  TCF,n  and 
TCF'2’  were  plotted.  Important  SAW  device  design  parameters,  such  as  coupling  coefficient,  velocity  and 
power  flow  angle,  were  computed  to  characterize  each  area.  As  part  of  Task  II,  experimental  results 
establishing  the  degree  of  correlation  with  theory  were  obtained. 

Theoretical  calculations  have  been  in  good  agreement  with  experimental  results.  Doubly  rotated  cuts 
of  quartz  with  an  improvement  of  TCF(2)  by  at  least  a  factor  of  two  over  the  ST  cut  have  been  obtained.  A 
further  improvement  was  obtained  after  a  second  iteration. 

Paragraph  2  presents  experimental  determination  of  propagation  characteristics  which  illustrate  an 
excellent  agreement  between  theoretical  calculation  and  experimental  results.  It  also  contains  the  results  of 
the  doubly  rotated  cut  TCF  measurements  made  to  date.  In  Paragraph  3  the  second  iteration  theoretical 
results  are  presented.  Doubly  rotated  cuts  with  superior  temperature  characteristics  to  those  discovered  in 
the  first  iteration  are  presented.  In  Paragraph  4  the  temperature  variation  of  the  power  flow  angle  and  a  new 
mask  design  insensitive  to  this  variation  are  discussed. 

2.  EXPERIMENTAL  RESULTS 

In  Section  III,  experimental  results  were  presented  demonstrating  the  good  agreement  with  the 
theory.  The  results  of  the  experimental  work  performed  are  summarized  in  Table  9.  The  orientations  of 
Table  9  are  expressed  in  the  notation  of  the  IRE  1949  standards*  and  were  obtained  by  selecting 
measurements  with  zero  first-order  temperature  coefficients  of  frequency.  This  was  done  to  ease  the 
comparison  between  the  experimental  and  theoretical  results  because  these  angles  can  be  accurately 


’"Standards  on  Piezoelectric  Crystals  1949”,  Pro.  IRE  14,  Dec.  1949,  pp.  1378-1395. 
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calculated  theoretically,  and  because  it  is  a  necessary  condition  for  low  TCF  orientation.  A  more  accurate 
measurement  of  the  second-order  temperature  coefficient  of  frequency  usually  results  when  the  first-order 
coefficient  is  small.  The  experimental  and  calculated  results,  as  can  be  seen  from  Table  9.  have  maintained 
their  high  degree  of  accuracy. 

Figure  35  shows  a  typical  frequency-temperature  curve  used  to  generate  the  data  in  Table  9  The 
stars  represent  experimental  data  points.  The  solid  lines  are  linearly  regressed  curves  used  to  define  the 
measured  first-,  second-,  and  third-order  TCF's  for  these  cuts  found  in  the  table.  The  first  order  temperature 
coefficient  is  small  enough  so  that  only  a  light  mask  rotation  is  required  to  arrive  at  an  orientation  with  a 
zero  TCF".  This  small  rotation  is  not  enough  to  significantly  change  the  second-order  TCF.  Experimental 
devices  with  a  near  zero  first-order  TCF  and  a  second-order  TCF  of  approximately  -1 .0  x  10S/C°2  have  been 
measured.  These  measured  second-order  TCF’s  are  significantly  lower  than  the  lowest  previous 
measurements  of  approximately  -1.5  x  10'8  with  a  zero  TCF1".  Figures  36  through  41  are  graphs  of 
experimental  frequency  temperature  characteristics  of  devices  fabricated  at  optimum  orientations. 
Changes  in  the  frequency  behavior  as  a  function  of  mask  alignment  can  be  seen  clearly. 

TABLE  9.  COMPARISON  OF  EXPERIMENTAL  AND  CALCULATED  RESULTS 


ANGLES 

CALCULATED 

MEASURED 

PHI 

THET 

E9 

TCF’tf 

TCF'ttt 

TCF*ttt 

TCF*1tt 

TCF’ 

TCF* 

TCP* 

m 

43.75 

ot 

’0.07  X  10-* 

0.00  X  10-* 

-o.«o  x  io -T 

0.11  X  io-’# 

-0.1  x  io-‘ 

-0.J7  X  10-7 

-0.17  X  10r1® 

K 

25.9 

135.7 

-0.01 

0.74 

-0.15 

0.42 

0.16 

•0.16 

0.54 

26.66 

1)4.9 

-0.24 

0.55 

-0.15 

0.43 

0.026 

-0.16 

0.47 

o 

26.97 

134.6 

-0.10 

0.60 

-0,15 

0.46 

0.067 

-0.13 

0.46 

MSI 

27.63 

134.2 

-0.26 

0.54 

-0.15 

0.49 

-0.06 

-0.15 

0.63 

6  00 

26.97 

135.6 

-0.04 

0.75 

-0.14 

0.46 

0.15 

-0.13 

0.30 

s.so 

27.63 

135.2 

-0.15 

0.65 

-0.14 

0.49 

0.13 

-0.13 

0.26 

-0.03 

26.70 

136.9 

0.29 

1.10 

-0.11 

0.44 

0.43 

-0.11 

0.49 

-  0.13 

29.1 

1377 

0.06 

0.69 

-0.11 

0.53 

0.24 

-0.097 

0.25 

0.633 

26.15 

137.0 

-0.24 

0.57 

-0.13 

0.47 

-0.11 

-0.15 

0.06 

-0  967 

26.23 

136.4 

0.01 

0.62 

-0.12 

0J1 

0.34 

-0.13 

0.79 

tWAFER  OBTAINED  COMMERCIALLY.  ANGULAR  TOLERANCE  IS  UNKNOWN. 
ttCALCULATED  USING  SINHA  AND  TIERSTEN'S  RROGRAM. 
tttCALCULATED  USING  FINITE  DIFFERENCE  ARRROACH. 

3.  THEORETICAL  RESULTS 

The  first  iteration  results  were  utilized  to  refine  the  theoretical  results  obtained  earlier.  The  close 
agreement  between  theory  and  experiment  demonstrated  earlier  prompted  a  more  thorough  theoretical 
search  for  promising  temperature-stable  doubly  rotated  orientations.  An  attempt  to  determine  optimum 
orientations  was  felt  justified  in  light  of  the  accuracy  of  the  calculated  results  Second-order  TCF's  :re 
predicted  by  the  finite  difference  method  to  within  0.005  ppm/COJ.  First-order  TCF’s  are  predicted  to  within  3 
ppm/C°  by  the  Sinha  and  Tiersten  perturbation  method. 
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Figure  35.  Frequency-Temperature  Dependence  for  (YX  wit)  0.633/26.15/137.0 

A  more  thorough  search  in  and  near  previously  discovered  near-optimum  orientations  did  not  yield 
any  new  results.  Other  areas  previously  considered  as  not  quite  optimum,  on  closer  examination,  were 
found  to  have  a  greater  predicted  temperature  stability  than  any  other  orientations  yet  measured.  A  second- 
order  TCF  of  -0.93  x  10®/C°2  with  a  zero  first-order  TCF  is  predicted  for  a  new  family  of  cuts.  The 
significance  of  this  family  of  temperature-stable  cuts  will  become  evident  only  after  devices  are  built  and 
experimental  results  are  compared,  as  the  difference  in  TCF,2,'s  is  about  10  percent.  The  addition  of  this 
family  of  cuts  to  our  investigation  is  also  significant  as  it  opens  up  the  possibility  of  selecting  cuts  not  only 
for  improved  temperature  stability  but  for  improved  stress  compensation. 

Table  10  contains  a  list  of  crystal  orientations  with  zero  first-order  TCF’s  and  low  second-order  TCF’s. 
Two  new  families  of  cuts,  centered  about  (YX  wit)  15o/30°/38°  and  (YX  wit)  12.5°/35°/130°  are  included  in 
Table  10.  Second-order  TCF,2,,s  of  -1.0  x  10"®  and  -0.93  x  10"®  are  predicted  for  the  two  families.  Each  of  the 
families  extends  over  a  surface  with  a  PHI  variation  of  20°  and  15°  respectively,  and  a  THETA  variation  of 
about  5°  and  10°  respectively,  the  PSI  angle  for  obtaining  a  zero  TCF'2’  being  fixed  for  each  PHI  and  THETA. 
Figures  42  through  45  show  these  areas  in  detail.  Both  the  zero  TCF"’  curves  calculated  by  Sinha  and 
Tiersten’s  perturbation  approach  and  the  finite  difference  method  are  shown.  Overlapping  second-order  TCF 
contours  are  also  plotted. 
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Figure  37.  Frequency-Temperature  Dependence  for 

(YX  wit)  0.633/26.15/136.813  to  (YX  wit)  0.633/26.15/137.016 


TEMPERATURE  °C 

Figure  38.  Frequency-Temperature  Dependence  for 

(YX  wit)  5.583/27.833/134.994  to  (YX  wit)  5.583/27.833/ 135.794 


TEMPERATURE  °C 

Figure  40.  Frequency-Temperature  Dependence  for 

(YX  wit)  7.417/27.833/134  to  (YX  wit)  7.417/27.833/135.2 
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TEMPERATURE  °C 

Figure  41.  Frequency-Temperature  Dependence  for 

(YX  wit)  8.033/26.9667/134.618  to  (YX  wit)  8.033/26.9667/136.818 


TABLE  10.  PROPAGATION  CHARACTERISTICS  OF  CRYSTAL  ORIENTATIONS 


0.4  X  10-® 
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Figme  45.  TCP  Contour  Map  (PHI  =  30°) 


Table  11  contains  some  of  the  propagation  characteristics  of  these  selected  orientations.  The 
velocities,  coupling  coefficients,  and  power  flow  angles  are  seen  to  be  quite  reasonable.  Tables  12  and  13 
contain  the  derivatives  of  the  first-order  temperature  coefficients  of  frequency  with  respect  to  angular 
changes  in  orientation.  These  derivatives  determine  the  accuracy  with  which  one  must  fabricate  to  ensure  a 
negligable  TCF"1.  The  nonzero  value  of  9TCF("/9<*>  allows  a  compensation  for  misalignment  of  the  wafer 
orientation  by  a  compensation  in  mask  orientation  as  the  values  of  9TCF 9TCF wldO,  and  9TCF,2l/9^ 
are  small.  In  the  first  iteration  of  Task  II,  it  was  demonstrated  that  devices  with  these  TCF(1>  angular 
derivatives  are  easily  fabricated  with  low  TCF1".  These  two  families  of  cuts  are  believed  to  be  entirely 
suitable  to  temperature-stable  SAW  device  fabrication. 

TABLE  11.  PROPAGATION  CHARACTERISTICS  OF  SELECTED  ORIENTATIONS 


MMLit  of  rrcfOi,  Momu 

It  AMO  Tt  fWOOWAM) 

VELOCITY 

K» 

POWER  PLOW 

TNVTA 

PH 

• 

2S 

128.21 

3288.84 

1.12 

-0.3 

• 

27 

12S.82 

3283.80 

1.12 

-0.2 

• 

2S 

128.58 

3280.83 

1.12 

-0.1 

7 

24 

126.88 

2203.32 

1.12 

-0.5 

7 

27 

12S.84 

3288.70 

1.12 

-0.4 

7 

24 

125.27 

3288.23 

1.12 

-0.3 

t 

2S 

125.74 

2310.15 

1.12 

-0.7 

• 

27 

121.1* 

3208.11 

1.12 

-0.8 

• 

24 

■m.tT 

3302.22 

1.10 

-0.5 

1 

24 

117.7* 

3288.80 

1.10 

♦0.7 

1 

27 

117.4* 

12*7.44 

1.10 

♦0.8 

1 

24 

117.17 

2288.28 

1.10 

♦1.0 

0 

24 

13*07 

3284.08 

1.12 

♦0.8 

0 

27 

117.7* 

3283.08 

1.10 

♦1.1 

0 

28 

117.** 

3282.28 

1.10 

♦1.2 

-1 

2S 

128.27 

2288.85 

1.10 

♦1.1 

-1 

27 

128.08 

3288.01 

1.10 

-1.3 

-1 

2S 

127.80 

3288.84 

1.08 

♦1.5 

14 

28 

40.188 

2288.80 

0.88 

-7.7 

14 

40 

40.418 

3208.87 

0.88 

-8.1 

14 

41 

40.84 

3318.18 

0.84 

-8.8 

IS 

28 

1*.7* 

3301.82 

0.88 

-78 

IS 

40 

40.00 

U10.U 

0.M 

-88 

IS 

41 

40.22 

3318.08 

0.88 

-8.8 

IS 

28 

28.4 

3306.28 

0.88 

-8.0 

IS 

40 

KMC 

3314.03 

0.88 

-8.4 

7.S 

28.0 

41.77 

3282.43 

1.00 

-4.8 

IS 

2S.0 

40.82 

13*4.21 

1.01 

-6.3 

128 

20.0 

M.M 

124127 

1.01 

-4.1 

12.S 

328 

28.4 

3284.48 

1.00 

-4.8 

1S.S 

20.0 

28.12 

3244.88 

0.88 

-4.8 

1S.0 

328 

28.55 

3287.82 

0.08 

-68 

TABLE  12.  aTCF"79*  FOR  SELECTED  CUTS 


[£ 


s 


Angles  of  ZTCFm,  Degrees 

(S  and  T'S  Progran 

1) 

aTCF'/a* 

PHI 

THETA 

PSI 

15.0 

32.5 

38.55 

-3.6 

17.5 

30.0 

37.35 

-3.4 

20.0 

30.0 

36.6 

-3.5 

22.5 

30.0 

35.85 

-3.5 

25.0 

30.0 

35.07 

-3.3 

27.5 

30.0 

34.28 

-3.4 

7.5 

35.0 

132.68 

*2-7 

10.0 

40.0 

129.40 

+2.8 

12.5 

35.0 

130.62 

+2.0 

425 

45.0 

127.15 

124.43 

+2.5 

+2.6 

17.5 

40.0 

126.14 

+2.1 

20.0 

40.0 

124.92 

+2.6 

TABLE  13.  aTCF("/a^  AND  aTCF"’/ae  FOR  SELECTED  CUTS 


Angles  of  ZTCF"’,  Degrees 
(S  and  T's  Program) 

aTC F"Va* 

aTCF'Vae 

PHI 

THETA 

PSI 

7 

27 

135.64 

-0.7  (PPM/C°)/DEGREE 

-0.5  (PPM  /  C°)  /  DEGREE 

0 

27 

137.78 

-0.8 

-0.8 

15 

40 

40.00 

+1.5 

-0.7 

15.0 

32.5 

38.55 

-1.2 

0.6 

12.5 

35.0 

130.62 

0.95 

_ i 

1.1 

4.  TEMPERATURE  VARIATION  OF  THE  SAW  POWER  FLOW  AN6LE 

An  unexpected  rapid  fluctuation  of  the  surface  acoustic  wave  power  flow  angle  on  doubly  rotated  cut 
quartz  saw  devices  was  discovered  during  the  testing  of  temperature-stable  SAW  devices.  In  this  paragraph 
the  phenomenon  of  a  temperature  variation  in  the  SAW  power  flow  angle  will  be  discussed. 

The  power  flow  angle  for  a  particular  direction  of  propagation  is  an  important  parameter.  While  the 
phase  fronts ‘always  remain  parallel  to  the  source  transducer,  the  wave  as  a  whole  does  not  propagate 


147 


perpendicular  to  the  wavefronts.  This  is  a  characteristic  of  anisotropic  substrates  where  the  phase  velocity 
is  asymmetric  about  the  propagation  direction;  i.e.H^ + Ai/»)  *  i #  -  W).  The  major  problem  which  arises  is 
that  the  acoustic  beam  may  steer  off  the  desired  propagation  track  missing  the  output  transducer  unless  it  is 
properly  designed. 


The  power  per  unit  width  carried  in  a  surface  wave  is  found  by  integrating  the  mechanical  and  electrical 
Poynting  vectors,  to  obtain 


where  m  is  the  particle  displacement,  T,j  the  stress  tensor,  <f>  the  electric  potential  and  D,  the  electric 
displacement.  Pi  and  P2  give  the  power  flow  perpendicular  and  parallel  to  the  wavefront,  respectively.  P3 = 0 
for  the  Rayleigh  wave  which  is  confined  to  the  surface.  The  power  flow  angle  may  be  defined  as  e  =  arctan 
(P2/P1).  Power  flow  angles  as  high  as  20°  are  not  uncommon  on  quartz. 


Figure  46.  Power  Flow  Angle  Against  Temperature  for  (YX  wit)  14  283/39-117/40-6 
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Calculations1*4  of  the  power  flow  angle  at  different  temperatures  for  doubly  rotated  cuts  of  quartz 
have  been  performed.  The  cuts  are  designated  by  the  1949  IRE1  standard.  Figure  46  illustrates  the 
temperature  variation  of  the  power  flow  angle  for  the  doubly  rotated  cut  of  quartz  (YX  wit) 
14-283/39-1 17/40-6.  The  important  feature  of  this  dependence  is  the  large  variation  of  the  power  flow  angle 
over  the  temperature  range  shown. 

Figure  47  contains  a  pictorial  representation  of  a  device  fabricated  at  (YX  wit)  14  3/39-1/40-6.  The 
input  transducer  on  the  left  generates  an  acoustic  wave  which  only  partially  illuminates  the  output 
transducer  on  the  right.  Figure  48  (a  and  b)  are  photographs  of  the  device  response  with  a  short  gated  RF 
pulse  as  the  input,  showing  the  response  at  131°C  and  34°C,  respectively.  The  first  notch  is  a  result  of 
missing  finger  pairs.  The  anisotropy  parameter5  was  calculated  to  be  0-625  at  -50°C,  0-614  at  25°C  and  0-586 
at  150°C.  The  transducer  apertures  are  34  mils  and  24  mils,  the  length  of  the  device  is  260  mils  and  the 
acoustic  wavelength  is  0-48  mils.  The  temperature  dependent  effects  of  diffraction  on  the  envelope  were 
found  to  be  negligible.  The  shortening  of  the  device  response  is  clearly  evident  from  the  photographs  and  is 
due  to  the  rapidly  decreasing  power  flow  angle  successively  illuminating  more  of  the  output  transducer  as 
the  temperature  increases.  The  power  flow  angles  estimated  from  these  photographs  are  displayed  in  Figure 
42,  alongside  the  theoretical  results.  Table  14  contains  a  summary  of  the  power  flow  angles  temperature 
dependence  for  several  temperature-stable  doubly  rotated  cuts. 


B.  IMPULSE  RESPONSE  AT  150°C 

C .  IMPULSE  RESONSE  AT  25°C  4453  ,3 

0.  IMPULSE  RESPONSE  AT  -50°C 

/ 

Figure  47.  Pictorial  Representation  of  Device  Response 

"‘Higher  Order  Temperature  Coefficients  of  the  Elastic  Stiffness  and  Compliances  of  Alpha-Quartz", 
Bechmann,  R.,  Ballato,  A.,  and  Lukaszek,  T„  IRE  Trans.,  1962,  pp.  1812-1822. 

’“Compensation  for  Oiffraction  in  SAW  Filters",  Savage.  E.  B.,  and  Matthaei,  G.  L,  1979  IEEE  Ultrasonics 
Symposium,  CH1482-9/79/,  pp.  527-532. 
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When  the  oscillator  frequency  is  not  exactly  at  the  synchronous  frequency  of  the  SAW  del:-,  h.  .0  the. 
phase  response  of  the  delay  line  becomes  temperature  dependent.  The  use  of  doubly  rotated  cut  SAW 
devices  clearly  requires  designs  which  can  accommodate  a  large  variation  in  the  SAW  power  flow  angle 
Consequently,  a  special  mask  used  for  testing  was  designed. 

Figure  49  contains  an  illustration  of  the  pattern  to  be  used  for  doubly  rotated  cuts  with  a  large  power 
flow  angle  temperature  variation.  The  large  aperture  of  the  output  transducer  is  required  to  accommodate 
rapid  variation  of  the  power  flow  angle. 

The  center  frequency  of  the  device  is  approximately  260  MHz.  The  electrode  and  line  spacing  are  0.06 
mil  (double  electrodes),  so  the  acoustic  wavelength  is  0.48  mil.  The  input  transducer  has  40  pairs  of 
electrodes  and  the  output  transducer  has  11  sections  of  electrodes  with  4.5  pairs  in  each  section.  The 
spacing  of  each  section  is  19.2  mils.  The  width  of  the  input  transducer  is  36  mils,  and  that  for  the  output 
transducer  is  120  mils.  The  design  allows  a  maximum  of  12  degrees  beam  steering,  which  is  adequate  for 
most  of  the  desired  cuts  in  this  study.  The  final  mask  is  a  stepped  and  taduced  design  having  a  family  of 
patterns  with  relative  rotation  of  a  fraction  of  a  degree,  so  that  a  family  of  TCF’s  with  small  increments  of  the 
propagation  angle  0.4°  can  be  measured. 

The  input  transducer  is  divided  into  two  sections.  Selecting  the  upper  or  lower  section  ensures  that 
the  complete  wavefront  of  the  acoustic  waves  will  propagate  over  the  output  transducer  aperture 
throughout  the  temperature  range  of  interest. 

In  designing  filters  and  reflective  array  devices  on  rotated  cuts  of  quartz,  additional  care  must  be 
taken  to  ensure  that  the  temperature  dependent  power  flow  angle  does  not  degrade  device  response.  If 
improperly  designed,  device  time  delay,  bandwidth  and  phase  all  become  temperature  dependent  when 
fabricated  on  a  rotated  cut.  These  observations  are  especially  important  for  reflective  array  devices  in 
which  rotated  cuts  of  quartz  are  often  used  to  achieve  temperature  compensation  in  two  different 
propagation  directions.  Suitable  device  design  can  overcome  these  problems  and  result  in  temperature 
compensated  delay  lines  useful  for  oscillator  applications.  The  moving  acoustic  beam  may  be  made  to 
illuminate  different  parts  of  the  output  transducer  at  different  temperatures  in  such  a  way  as  to  maintain  a 
steady  phase  over  the  entire  temperature  range.  A  design  using  the  temperature  dependent  PFA  is  currently 
being  fabricated  in  which  a  compensation  signal  maintains  a  relatively  constant  output  phase  despite  the 
temperature  dependent  delay  time.  Input  beam  aperture  compression  may  be  used  to  achieve  increasingly 
accurate  adjustments  of  the  output  phase.  Materials  such  as  lithium  niobate  with  higher  coupling 
coefficients  may  be  useful  in  this  respect. 


Figure  49.  Transducer  Design 
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By  proper  selection  of  the  output  transducer  location,  one  may  obtain  temperature  stable  SAW 
oscillator  due  to  the  power  flow  angle  sweeping  as  a  function  of  temperature.  Figure  50  shows 
compensated  SAW  oscillator  using  the  combination  of  two  90°  out  of  phase  transducers  shown  in  Figure  49. 
Where  the  combined  output  phase  of  the  transducers  is  compensated  to  a  minimum  change  as  the 
temperature  varies;  hence  improved  temperature  characteristics  are  obtained. 

5.  REPRESENTATIVE  EXPERIMENTAL  DATA  IN  THE  SELECTED  AREAS 

Several  selected  areas  have  zero  first  order  TCF  and  low  TCF'2’  and  TCF(3)  have  been  fabricated  and 
evaluated.  These  areas  provide  some  of  the  best  results  in  the  study.  The  cut  YX  wit  <*>/0/ ifr  of 
12.01/31.21/39.8  has  less  than  80  PPM  shift  of  frequency  in  the  temperature  -40°C  to  +140°C.  The 
experimental  result  is  shown  in  Figure  51 .  This  cut  has  measured  TCF"’  of  0.55  x  10  7  TCF'21  of  -1 .19  x  10  s  and 
TCF131  of  0.45  x  10  10  Figure  52  shows  the  experimental  results  of  a  device  in  the  15.13/31.13/37.33 
orientation.  The  device  has  TCF"’  of  1 .1  x  10'7,  TCF(2)  of  -1 .18  x  10~8  and  TCF131  of  0.253  x  10~‘°.  This  device  has 
less  than  50  PPM  shift  from  -30°C  to  +100°C  temperature  range.  Other  selected  cuts  are  shown  in  Figures  53 
and  54.  These  cuts  all  have  relatively  good  temperature  characteristics. 
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Figure  50.  Power  Flow  Angle  Compensated  SAW  Device 


Figure  52.  Temperature  Stable  Cut  of  Quartz  for  SAW 
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Figure  53.  Temperature  Stable  Cut  of  Quartz  for  SAW 


SECTION  V 
CONCLUSION 


The  temperature  coefficients  of  frequency  (TCF)  have  been  analyzed  for  doubly  rotated  cuts  of 
quartz  for  surface  acoustic  wave  devices.  The  analysis  procedure  is  shown  below: 

(1)  Using  Finite  Difference  program,  identify  angular  areas  with  zc  _  P"  and  low  TCF12’  on 
10°  x  10°  x  10°  grid. 

(2)  Calculate  in  these  areas  with  2.5°  x  2.5°  x  2.5°  grid  to  locate  the  rotation  with  minimum 
TCF'2’  which  has  TCF"’  =  0. 

(3)  Identify  the  accurate  zero  TCF11'  with  Sinha  and  Tiersten's  program  in  areas  obtained 

in  (2). 

(4)  Calculate  TCF'1',  TCP2’  and  TCP”  in  selected  areas  obtained  in  (3)  with  1°  x  1°  x  1°  span. 

(5)  Calculate  the  coupling  coefficients,  propagation  velocities,  power  flow  angles  and 
inverse  bulk  wave  velocities  for  these  angles. 

The  conclusions  of  the  study  are: 

(a)  The  zero  TCF"1  surface  and  zero  TCP2’  surface  data  does  not  show  an  intersection 
in  the  10°  x  10°  x  10°  grid.  However,  orientations  with  frequency  variation  <*>  100  ppm  from  -50°C  to  100°C, 
better  than  ST  Cut  by  a  factor  of  two ,  have  been  found  in  several  areas  where  TCP 1  ’ = 0  and  TCP21  is  less  than 
-1.5  x  10=«/Co  J. 

(b)  The  aTCF/a  angle  is  normally  very  large  in  the  areas  with  zero  TCP"  and  low 
TCP2',  therefore  accurate  crystal  orientation  is  critical.  It  is  estimated  that  the  orientation  accuracy  of  6 
minutes  is  required  to  obtain  a  TCP"  of «  50  ppm  in  the  -50°C  to  100°C  temperature  range  in  the  three  areas 
investigated.  These  are  <t>lol\b of  6/27/135.93  (4.45  x  10'2  ppm/°C  •  minute),  14/40/40.415  (5.83  „  10~2  ppm/°C 
•  minute)  and  0/27/137.78  (5.00  x  10~2  ppm/C°  ■  minute).  Orientations  with  less  sensitivity  to  orientation 
accuracy  have  been  evaluated.  These  are  the  cuts  that  have  3[TCP"]/di)r  =  0  when  TCP"  =  0.  The  cuts 
evaluated  were  =  20/30/155,  20/20/150,  10/40/168.  The  TCP2’  of  these  cuts  are  approximately 
4  x  10 "/C0,  comparable  or  higher  than  ST  cut. 

(c)  Coupling  coefficients,  velocities  and  beam  steering  angles  were  defined  for  the 
selected  areas,  no  leaky  modes  were  found  in  these  areas.  Determinations  of  experimental  propagation 
characteristics  using  the  laser  probe  technique  demonstrated  good  agreement  between  measured  and 
calculated  values. 
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(d)  Methods  to  X-ray  orient  the  quartz  crystals,  cut  the  doubly  rotated  wafers  and 
polish  such  crystal  wafers  were  developed. 

(e)  It  has  been  determined  that  TCP"  is  the  most  sensitive  to  angular  variations; 
thus  the  large  values  of  dTCFI/^  imposes  strict  fabrication  tolerances  on  processing.  It  was  necessary  to 
incorporate  a  mask  design  using  rotated  structures  to  compensate  for  processing  errors  within  the  allowed 
fabrication  tolerances. 

(f)  Theoretical  c  ilculations  of  power  flow  angle,  coupling  coefficients,  and  second 
and  third  order  TCP’s  indicate  that  they  do  not  vary  quickly  with  angle.  Thus,  by  varying  <h  on  any  particular 
cut,  linear  temperature  variations  can  be  compensated  for  while  all  other  cut  parameters  vary  only  slightly. 

(g)  Experimental  design  and  procedures  followed  during  the  measurements  of  the 
TCFs  indicate  that  the  experimental  error  is  within  10  PPM.  Tables  8  and  9  which  display  the  comparison  of 
theoretical  and  experimental  results,  illustrate  the  excellent  agreement. 

(h)  Experimental  results  of  the  1st  iteration  indicate  that  an  improvement  in  TCP" 
by  at  least  a  factor  of  two  over  ST  quartz  can  be  obtained.  A  second  iteration  approach  using  a  higher 
resolution  and  better  selected  areas  has  improved  the  TCP"  of  the  doubly  rotated  quartz  over  the  ST  cut  by  a 
factor  of  three  or  better. 

The  best  experimental  data  was  obtained  with  the  devices  in  the 
12.01/31.21/39.8  orientation.  The  total  shift  of  frequency  of  less  than  80  PPM  in  the  -45°C  to  +140°C  was 
obtained.  Several  other  areas  that  will  provide  comparable  temperature  characteristics  are  described. 

(i)  T"°  temperature  variation  of  the  power  flow  angle  on  doubly  rotated  cuts  of 
quartz  and  its  effect  on  device  design  was  discussed.  An  oscillator  design  accommodating  large 
temperature-varying  power  flow  angles  was  discussed.  The  utilization  of  the  temperature  variation  of  the 
power  flow  angle  to  achieve  temperature  compensation  of  the  oscillator  was  demonstrated. 
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APPENDIX  A 

VOLUME  PERTURBATION  OF  AULD 


Perturbation  techniques,  as  used  here,  allow  calculations  of  small  changes  in  the  solutions  to  a 
problem  caused  by  small  changes  in  the  physical  parameters  of  the  problem,  once  the  solution  to  the 
unperturbed  problem  is  known. 

Our  general  approach  to  the  problem  of  determining  the  temperature  dependence  of  V.  will  be  as 
follows.  First,  the  Rayleigh  wave  propagation  problem  will  be  solved  in  the  standard  way  in  its  entirety  at 
room  temperature,  T0.  Given  the  solution  of  problem  at  T0  and  the  dependence  of  the  physical  constants 
(such  as  c„)  on  temperature  at  T0,  we  will  apply  the  volume  perturbation  formula,  calculating  the 
temperature  dependence  of  Vs.  The  dependence  of  V,  on  T  is  then  used  to  calculate  the  frequency 
characteristics  of  the  actual  device  given  the  thermal  expansion  coefficients  as  a  function  of  temperature.  At 
this  point,  the  frequency  temperature  dependence  of  the  substrate  as  a  function  of  crystal  cut  and  direction 
can  be  thoroughly  explored. 

The  Volume  Perturbed  Formula. 

Denoting  u,  (T0)  by  u„  etc.,  the  volume  perturbation  formula  is  given  by1: 

=  [A]  dy if  (~u*T  -u'T* 

O  O 

+  0*  (icuD')  +  (i<uD)*)  •  z  dy  (A-1) 

[A]  =  [Ap  u*  •  u'  +  T*:  (ase:  T  +  Ad  ■  E') 

+  E*  •  (A*t  •  E'  +  Ad:T')] 

This  equation  is  exact,  but  involves  knowing  the  solution  u'.  However,  we  can  set  u  =  u'  if  the  temperature 
dependence  of  the  solution  is  small  and  use  P  =  (1/2)  Re  J  (-u*T  +  E  H*)  dy  to  obtain  the  approximate 
solution 

A0r  =  (w/4P)  J*  [Ap  |  u  |  1  +  S:  Ac:S 

v  o  (A-Z) 

+  E*  •  a«t  •  E  +  E*  •  Ae:S  +  S:Ae  •  E]  dy 
The  u's  and  E's,  and  p  come  directly  from  the  computer  solution  at  T0. 

'Auld,  B.A.,  “Acoustic  Fields  and  Waves  in  Solids,"  Vol.  II. 
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The  a  terms  can  be  written  as  expansions  of  the  form: 


A0R  =  0*  (To Map'"  dT  +  a/'  dT2  +  aj"  dT3) 

Ap  =  p(T0)(a(p  dT  +  «<2>  6V  +  Q«;  dT3) 

=  p'"  dT  +  p'2'  dT2  +  p'3'  dT3 
ACiiki  =  c«ki'"  dT  +  c,)kl<2'  dT2  +  c,jki,J'  dT3,  etc. 

to  obtain 

Aj8„  =  j8r  (To)  ( afi'"  dT  +  /3r  (To)  6V  +  /SR  (T.)  a^'3'  dT3) 
=  MT/4P)  J*“  (p  (To)  |  u  | 2  +  S:c‘":S  +  E*  • •  E 

O 

+  E*  •  e"’  •  S  +  S:e‘"  •  E)  dy 
+  MT2/4P)  f  (p(T0)  a  <2'  |  u  |  2  +  S:cl2':S  + ...)  dy 

O  H 

+  MT74P)  f  (p(T0)  Op'31  |  u  !  2  +  S:c,5,:S  +  ...)  dy 

O 

or 

V  =  (<o/4P/3r  (To))  f  (p(To)  |  u  | 2  +  S:c“':S  + ...)  dy 


(A-3) 

(A-4) 

(A-5) 

(A-6) 


(A-7) 


(A-8) 


For  the  problem  of  quartz,  the  electrostatic  coupling  is  small  and  the  electric  terms  can  be  ignored. 
Thus  Ad  and  a<  in  (A-8)  may  be  set  to  zero  resulting  in  the  simplified  equation. 

=  (W4P/J*)  X°°  (p(T0)  «*  |  u  | 2  +  S :  cu>:  S)  dy  (A-9) 

o 

The  integrals  are  calculable,  for  instance,  from  (14), 
p  a->  i  o  i  *  a» 

oo 

=  p(T,)  «'p  5  /„  [£  C™  a,""'  exp(— i  0*  b'""  y)]  (A-10) 

[2  C;  a,'"'*  exp  (i0Rb""  y)*]  dy 
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The  problem  can  easily  be  extended  to  multi-layered  media  by  simply  performing  the  integral  over 
each  layer  separately.  For  the  double  layer  (A-9)  becomes: 

p,  (To)  <#;  |  u  | 1  dy  +  J*°°  P2  (T0)  |  u  | 2  dy  +  -  (A-11) 

o  a 

Because  the  dependence  of  u,  is  given  explicitly  by  (13),  it  is  not  necessary  to  use  numerical 
integration  procedures,  as  demonstrated  in  (A-10). 

Once  the  a ^  have  been  calculated  from  (A-9)  the  avu'  are  simply  determined.  Using  p  -  <u/V,,  and 
(A-9)  we  determine  aFU). 
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APPENDIX  B 


VOLUME  PERTURBATION  FORMULA  OF  SINHA  AND  TIERSTEN 

The  approach  used  by  Sinha  and  Tiersten  includes  the  effects  of  distortion  caused  by  heating  a 
substrate  of  quartz  and  is  based  in  the  coordinate  system  to  which  the  fundamental  elastic  constants  refer. 
In  this  reference  system,  the  density  is  constant  and  the  aF1'1  are  equivalent  to  the  av"1. 

In  this  new  reference  system,  the  perturbation  formula  becomes1 

(Av,/v.)  =  (1/2j8rV)(H/N2),  (B-1) 

4  4 

N2  =  (pVi/0*2)  2  X  (C(B'  a""1  Cw*  -  &*)  (B-2) 

m  =  I  n  =  1 

y.00  -  kI  fin 

H  =  -J  dy  J  dz  (K„  Ui.t  +  K„  uy„  +  Ku  u„  +  K„  u,.y  +  uy,y  +  K„  u„y)  (B-3) 

0  ~rrlfin 

KLy  =  (Cl7m«  "f"  ACL-yMa)  m  (B-4) 

For  the  first  order  perturbation  in  T  we  have: 

•*L>Ma  ~  (^l7mFab  °ab  ^l7km  a*>*.  +  ®yic)  ^T  (B-5) 

Ac(-yMfl(  =  (d  cLVM«/dT)  dT  v  (B-6) 

‘S 

Where  cLyM<*  are  the  second  order  elastic  constants  previously  denoted  simply  as  cLyM<*  and 
Ci  ywi-AB  are  the  third  order  elastic  constants.  The  terms  AcLyM<*  as  calculated  by  Tiersten2  is  available  only 

J  2 

to  first  order  in  T,  and  the  higher  order  elastic  constants  cLyMVABCD  have  never  been  determined. 

4 


"‘On  The  Temperature  Dependence  Of  The  Velocity  Of  Surface  Waves  On  Quartz,”  B.K.  Sinha  and  H.F. 
Tiersten,  1978  Ultrasoncs  Symposium  Proceeding,  IEEE,  pp.  662-665. 

"‘Temperature  dependence  of  the  Fundamental  Elastic  Constants  of  Quartz,"  B.K.  Sinha  and  H.F.  Tiersten, 
Proceedings  of  the  32nd  Annual  Symposium  on  Frequency  Control,  1978,  pp.  150-153. 
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APPENDIX  C 

THE  DIFFERENTIATION  METHOD 


A  method  for  determining  the  theoretical  temperature  dependence  of  Rayleigh  Surface  Waves 
consists  of  formally  differentiating  the  wave  equation  and  boundary  conditions.  The  boundary  conditions 
and  wave  equation  must  be  true  at  all  temperatures,  placing  restraints  on  how  the  parameters  of  the  wave 
equation  may  vary.  In  this  technique,  the  derivatives  of  these  equations  with  respect  to  temperature  are  set 
to  zero  and  solved  for  the  velocity  temperature  dependence.  This  method  follows  the  methods  used  by 
Bechmann,  Ballato,  and  Lukaszek1  to  compute  the  temperature  dependence  of  the  fundamental  elastic 
constants  from  frequency  data,  except  that  the  simplifying  assumptions  of  assuming  bulk  wave  solutions 
cannot  be  made.  This  method  was  later  used  by  Hauden2  to  search  for  temperature  stable  cuts  of  quartz. 

Christoffel’s  wave  equation  can  be  written  in  matrix  form  as 


*r„  -  pV.2 

r12 

r14" 

on 

rl2 

r22  -  pV.2 

r» 

r2« 

a2 

rl3 

r2J 

r»-pV,2 

rM 

ay 

r.4 

r2< 

Tl4 

r«4 

.a*. 

where 


=  0 


r„  =c5J  b2  +  2cub  +  Cn 

r2:  =  Cm  b2  +  2Cm  b  +  Cm 

rJ3  =  cm  b2  +  2c»  b  +  Css 
r  12  =  c<s  b2  +  (c  14  +  c$*)  b  +  Cm 

ru  —  Cjj  b2  +  (cu  +  Cjs)  b  +  Cu 


r23 

= 

Cm  b2 

+ 

(C36 

+ 

C45) 

b+cj« 

r44 

= 

~  («SJ 

b2 

+  2tu 

b  + 

til) 

ru 

= 

6js  b‘ 

+ 

(Bis 

+ 

831) 

b  + 

6|| 

r24 

= 

e»4  b2 

+ 

(6|4 

6m) 

b  + 

6i* 

rM 

= 

e»  b2 

+ 

(6ij 

+ 

63s) 

b  + 

615 

(C-D 


(C-2) 


'“Higher  Order  Temperature  Coefficients  of  the  Elastic  Stiffnesses  and  Compliances  of  Alpha  Quartz," 
Bechmann,  Ballato,  and  Lukaszek,  Proc.  IRE,  Aug  1962,  pp.  1812-1822. 

^‘Higher  Order  Temperature  Coefficients  of  Quartz  SAW  Ocsillators,”  D.  Handen,  M.  Michael,  J.J.  Gagnepain, 
Proc.  Frequency  Control  Symposium  (1978),  pp.  77-86. 


following  the  notation  used  in  the  previous  section  and  in  Matthews'.  This  equation  holds  for  jach  of  the  four 
modes,  hereafter  designated  by  a  superscript  m. 

The  boundary  conditions  become,  in  matrix  form, 


...(Cjjii  +  C3ii)  blm')  Oi1"’1  +  (6133  +  6333  blBI)  04<ml . 
...(C3111  +  C3113  b(ml)  «iln'>  +  (6131  +  6331  b1B))  o«ln” . 
{ Ci 2 1 1  +  C2313  b'ml)  +  (e,32  +  6332^"")  «4lm' . 

...(e„>  +  ew  b(m’)  oi ,m)-  (*„  +  *33  b,B)  -  i*„)  «4,B'. 


=  0 


(M) 


where  only  the  m’th  column  is  shown  and  the  Cm's  are  the  amplitudes  of  each  mode.  This  formulation 
assumes  a  normalization  of  the  a's 

(a,'"*)2  +  (02'"”)2  +  (a,""’)2  +  («4<m>)2  =  1.  (C-4) 


The  first  condition  we  can  place  on  the  wave  equation  is  that  the  determinant  of  Christoffel’s 
equation  vanishes  for  all  temperature,  or 


r..  -  pV,2 

rl2 

r.3 

r14 

r.2 

r2I  -  pV,2 

r2J 

r24 

r.3 

r23 

r„  -  pV,2 

r,4 

r,4 

r24 

r34 

r*4 

For  each  of  the  four  modes,  this  equation  is  valid,  resulting  in  four  equations  in  five  unknowns,  dVs/dT,  and 
db,m7dT  of  the  form  F,M  (dVs/dT,  db,m7dT)  =  0,  m  =  1,  2,  3,  4. 

For  each  mode,  Christoffel’s  matrix  equation  (C-1)  must  vanish,  yielding  the  set  of  equations 

(d/dT)  (a,"”'  (r„,m>  -  pV,2)  +  a2(n”  r,2,n”  +  as'"1  iV"'  +  a4,B’  r'w’)  =  0 
(d/dT)  (a,""'  r'u’+  a2,m'  (riT-  pV,2)  +  a,'B'  r,,'"’  +  a.(B>  rt? )  =  0 

(d/dT)  r„,n”  +  02,b>  r:3,B'  +  a,'-'  (r,3<Bl  -  pV.2)  +  rS? )  =  0  (C-6) 

(d/dT)  (a,,B’  r,4,n”  +  a2,m’  r24(B'  +  a3(0"  r.34<Bt  +  «4<B'  r44(B')  =  0. 

This  results  in  16  equations  (four  per  mode)  and  an  additional  16  unknowns  (the  oj|m,s)  of  the  form 
F2,'m>  (dV,/dT,  db(m|/dT,  d  ai,B,/dT,  d  «2,B'/dT,  d  a,'B)/dT,  d  a4,B7dT)  =  0 


"‘Surface  Wave  Filters,"  H.  Matthews,  John  Wiley  and  Sons,  New  York  (1977). 
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Using  (C-3)  we  obtain  the  single  equation 


(d/dt) 


(Cun  +  Cjji)b<m|)  (6uj  +  eujb'"1)  a/*1 
(Cjtii  +  Cji.jb'"”)  a,1"”  +  (6iji  +  ou1"’ 


=  0 


(07) 


of  the  form  F,|m'  (db(1|/dT,  db,27dT,..„  d«,,m7dT)  =  0 
and  four  equations  of  the  form 

(d/dT)  X  C™  [(Cjji,  +  c,wb'm»)  aT  +  (e.»  +  e„,b‘“')  aV"’  1  =  0 

(C-8) 

(d/dT)  X  Cm  [(cj,„  +  Cj,ijb<n") +  (e,„  +  eJ31b'-')  a/"'  ]  =  0 

m 

of  the  form  F^""'  (db(n”/dT,  da,'m7dT,  dCJdT)  =  0  (C-9) 

From  (5)  «,«■»  (d<*,""7dT)  +  a,'"1  (da^'/dT)  +  a,'"'  (da,'-7dT)  (C-10) 

+  a*"'  (da4m/dT)  =  0 
or 

Fj""1  (d  o,""7dT)  =  0  (C-11) 

Combining  the  above  results  gives  25  equations  in  the  25  unknowns. 
dV,/dT,  db""7dT,  dot'"7dT.  dC„/dT 

which  are  solved  simultaneously .  Once  the  temperature  dependence  ol  the  Rayleigh  wave  velocity  is  found, 
the  frequency  dependence  is  found  in  the  usual  way. 

To  obtain  the  second  order  dependence  of  V„  the  25  equations  are  differentiated  again.  The  values  of 
dV,/dT,  db(m7dT,  da,"”7dT  and  dCm/dT  previously  obtained  are  used  to  obtain  dJV,/dTJ,  etc. 
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APPENDIX  D 


SOLUTION  METHOD 

Calculating  the  velocity  of  a  Rayleigh  wave  requires  a  complete  solution  of  the  problem  to  be 
performed.  All  of  the  constants  in  equation  (13)  must  be  evaluated.  The  standard  used  at  Motorola  as 
developed  by  Campbell  and  Jones  is  outlined  below. 

First  the  fundamental  constants  are  rotated  into  the  coordinate  system  of  interest.  Next  a  value  of  V, 
is  picked.  The  coefficients  r,,  of  equation  (C-2)  are  evaluated  and  the  determinant  of  the  matrix  in  equation  (C- 
1)  is  set  to  zero,  as  it  must  be  if  a  solution  of  (C-1)  is  to  be  found.  This  results  in  an  eighth  order  equation  in  b. 
This  equation  is  solved  for  the  eight  complex  roots.  The  four  roots  b(n”  in  the  lower  complex  plane  are 
retained,  the  four  discarded  roots  not  satisfying  the  boundary  conditions  at  infinity.  Equation  (C-1)  is  then 
solved  for  the  four  eigenvectors  a,(m|.  The  b,"”>  and  the  corresponding  eigenvector  a,(m|  are  substituted  into 
equation  (C-3)  and  the  determinant  of  the  matrix  in  equation  (C-3)  is  evaluated.  This  determinant  must  be 
zero  for  a  solution  to  (C-3)  exist.  If  it  is  not  zero,  V.  is  varied,  and  the  whole  procedure  repeated,  until  it  is. 
Once  a  value  of  V,  is  found  such  that  the  determinant  in  equation  (C-3)  vanishes,  the  solution  to  (C-3)  is  found, 
giving  the  values  for  Cm.  These  constants  completely  describe  the  solution  of  equation  (13),  as  well  as 
providing  the  exact  velocity. 
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APPENDIX  E 


EQUIVALENCE  OF  TEMPERATURE  COEFFICIENTS  OF  FREQUENCY  AND  DELAY 

In  the  following  section,  the  relations  relating  the  first,  second,  and  third  order  temperature 
coefficients  of  delay  and  frequency  will  be  derived.  These  relations  show  implicitly  the  equivalence  of  each 
representation  of  the  device  temperature  characteristics,  and  justify  their  interchangeable  usage. 


In  the  text,  we  have  used  the  following  notation: 

t  =  delay  time  of  delay  line  oscillator  (E-1) 

ith  order  temperature  coefficient  of  delay  =  TCD"'  =  a"  (E-2) 

F  =  frequency  of  delay  line  oscillator  (E-3) 

ith  order  temperature  coefficient  of  frequency  =  of"’  (E-4) 

r  =  To  (1  +  a1;*  (T  -  To)  +  a,2;(T  -  T0)2  +  a  '3;.(T  -  T.) 5  +  ...)  (E-5) 

F  =  Fo  (1  +  of"’  (T  -  To)  +  oF'2’  (T  -  To)2  +  of'3’  (T  -  T.) 1  +  ...)  (E-6) 


where  T  is  temperature  and  T0  is  a  reference  temperature,  25°C  in.  our  case.  For  a  SAW  oscillator,  Ft  = 
constant  or 

rl  To  =  F„/F 

=  1/(1  +  of"’  (T  -  To)  +  of'2'  (T  -  To)2  +  of'5'  (T  •  To)3  +  ...)  <E'7) 


Using  the  relation 

1/(1  +  X)  a  1  -  X  +  X2  -  X3  for  X«1 


(E-8) 


we  can  write 


t/to  =  1  -  [of"'  (T  -  To)  +  of'2’  (T  -  To)2  +  oF'3’  (T  -  T.)3] 

+  [(of'"  (T  -  To)2  +  2of<"  (T  -  T.)  of'2’  (T  -  To)2] 

-  [of"1  (T  -  T„)]3  +  higher  order  terms  (E-9) 

=  1  -  of"’  (T  •  T.)  +  (-OF*2’  +  (of'")2)  (T  -  To)2  +  (-Of'3’  +  2of"’  of'2’  -  (op"’)3)  (T  -  T.)3  (E-10) 
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Equating  powers  of  T  -  To  in  (E-5)  and  (E-10),  we  obtain 


r"’=V" 

(E-11) 

r121  =W2’  +  (af"f 

(E-12) 

,;5'  =VJ,+2aFn'  of12’  - 

(E-13) 

As  the  only  assumption  on  F  and  r  used  is  that  of  equation  (E-7)  which  is  symmetric  in  F  and  r,  we 
immediately  obtain 


of*"-  -ar"’ 

(E-14) 

of'21  =  -a;2'  +  (ar">)2 

(E-15) 

of'31  =  W3’  +  2«rui  ar<2'  -  (arU1)3 

(E-16) 

Using  (E-11)  through  (E-16),  we  can  always  relate  one  set  of  temperature  coefficients  to  the  other. 
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APPENDIX  F 
X-RAY  RUN 

The  following  is  an  example  of  the  computer  results  from  the  FORTRAN  X-ray  orientation  computer 
program  for  the  doubly  rotated  cut  7/27/0.  File  MILDAT  contains  only  the  two  Miller  indices  shown  on  the 
printout  on  this  run. 
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APPENDIX  6 


X-RAY  PROGRAM 

The  following  is  the  FORTRAN  computer  program  used  to  calculate  X-ray  diffraction  angles  for  a 
doubly  rotated  cut  of  quartz. 


179 


pftBCEDIJtf)  PAflS 


u 

%u 

OL 

V 

2 


Z 

h  I 


2  UJ 

w  « 

2  iu 
2  ►- 
C  2 


xuz<< 

NXOWI 
— »  « 

z  z  m  o  i 
c  c  z 

Z  IhiJ 
S  Z  X  J 
o  O  M  «  < 
U  UO  U  l 


z 

c 

u 


C  I 
2  I 
w  • 


Ul  M 

#  2  e 

2 

M  W 

.  K  C 
I  2  < 
C  4J 

I  U  ft 


«  w  c 

J  • 
m  m  a 

X  K  W 
<K  O  • 
otft  c 
zra  m 
ui 

c 

«<kC 

u  u*-*o 


U  M  I 

©  <  : 


O 

o 

o 

■ 

a 

z 

Ui 


MO  /'M 

• Q  (5HK 
9  0  MM© 

C  ft  «  ^  vz 
*»lLt  h  K  * 


%  0)  0) 

tM  WM  M  «  * 

<UUU  MIW 
Ui  %*«%«)•) 

►>  •*  •-  MM  4  4 

4  «  «  0« 

►  \  •  %  •MO 

C  a:  ui  —  m  m  l  ♦  na 

a  N  J  4  4  4  MMMOC\| 

v  wwuuO  *  U  4  4  O 
4-  M  CZ  2  2  2  4  4 
2  «a4333— •Ai«**0U0O 
M  K>N  \  U  Ufe.  U0U  tO  f  •  I  U 
►-  w  uuc  «■■■■•>■• 


DU 

•  4 

OHV 
MM 
4  4 

MO 

•  ♦ 
MM 
MM 
4  4 


Z  J 

C  9 

o  u 


9-  C  , 
2  <  . 
C  Ui  • 
o  a  i 


l  C  2 

i »-  a 
C  c 
I  ©  u 


z  & 

2  O 

O  0- 

4.  00 


c 

a 

O  OD 

x: 

Ui  M 


-i  r 
0  2 

Ibj  c 

a  u 


-l  J 
J  -J 


-M  *»«Vi 

IUIU  fW  MM  < 

W  W  W  W  W  W  W  W  1 

laxtKtxa: 


ocooooooeoeoeocoeoccoceoeeoococeoc  c  o  o  o  o  0001 
occeeceoeeeeeeceeQeoeeeecoeeeooccc  o  o  c  o  ©  oo  e  i 
o  o  s  co  o  o  o  o ocmooccoooooco  oeoeooeooooooe  ©  c  e  ©  ©  i 


«.  «.  Ik  c  «.  •>  t 


i  ».  i  *.  *.  % 


4  I.  «  I.  c  I.  4  I. 


c  «•*<**>  c  —  aj  Kt «  an  «9  o-m^  « >  o 


*1.000  IfRR*! 

na'.OOO  UPTURN 


tel  • 

•  II 

*  te. 

%  » 

ft  -! 

«  % 

9-  tel 

X 

•5 

tef  ■ 

tel  •- 

te» 

X  *- 

*>  —l 

ft 

X  X 

•  X  ► 

► 

O  J 

V  C  *0 

% 

M0  OO 

•  ft  V 

X 

3  te. 

•  te  a 

0k 

X 

«  u 

UJ 

•  tel 

•M 

00  X 

ft  -l  O 

% 

13  •“ 

fO»- 

A* 

Z 

«  7 

w 

•m2 

«teK 

X 

•ft  O 

«  «0 

tel 

M 

ft  »  *M 

a 

tel  tel 

ft  o 

% 

X 

o 

ft  tel 

01  «A 

IA 

• 

ft  X  • 

M 

►  oft 

ft 

-1 

«  ft-  — 

% 

m 

ai 


v. 

•  z 
o  or 

•  z> 


ai 

a 


V. 

c 


x 

0) 


c 

o 


«  Jft. 
ft  ft 

•  X 

•  o-  o 

««K 

•  Ik 
«•- 

•  Ui 
*1  « 
«u« 

•  =>  * 
«•>  V 


tel 

te  X  mK 

•5 

ft  o 

•  M 

« 

a 

>  ft  tel 

•  ft 

ft  o  o 

ft  ft 

% 

ftM* 

OO  M 

► 

«  tel  -l 

J  tel 

M0 

off 

It*  O  M  MX 

9  M  M» 

< 

ft  2  U 

ft  to 

X 

O 

►  -I 

X  Al  • 

•  — 

ft 

•  •ft 

*  X 

10 

tei  * 

►  tel  ^  © 

ft  o  W 

ft 

«  te.  2 

MX  MX 

w 

Al 

x  te 

•>2  0 

te  K2 

•te 

ft  tel  o 

M  0ft 

X 

« 

*-  •» 

tel  3  ft  • 

►  C  ft 

% 

►*  ft  O  9te 

x  X 

tel 

•3  A> 

•  X 

u  ■  ► 

tel  CS  X 

ft 

ft  ft  ►“ 

—  Al 

ft 

H  2  m 

Z  V 

M>  VO 

tt  MX  MX 

O 

MO  ««3  ft 

ww 

3  V  X 

ft  ft 

m  m  m  • 

m  oA 

► 

ft 

O  ftklK 

X  X 

U  ft  ** 

a  z  tei 

>*  tel 

•-  o  ~ 

OO  W 

03 

X  ft 

MX  JO 

tel  tel 

tei  tel 

o  <  * 

X  -1 

►  Ute 

•  oZ 

MX 

ft.  ► 

O  OB 

•E  X 

too 

te.’  *m 

a  -i 

tel  2  tel  X 

»-  es  ft 

© 

X  a  -l-i  tel 

9  Z 

2  w  w  u  ft  •  A 

tel  -4 

C  *M 

m  i93  «  a 

-j  tei  — 

IT 

ft  ■  ■  9  ft 

9- 

•  ft  ft 

«Z  X 

JK  h  • 

2UZX 

2  3 

ADO  O  J 

•  ■  2 

Al 

R 

»OC  M 

**33 

ft  J  JO 

o 

ft 
c  « 


Z  AWZ 

O  2  13  C. 
X  Ui  to  3 

I  f »-  7 
OHZC 
UC  HU 


z  a 

c 

I  J 
7  ft 

C  U. 
u  * 


Al  W 

•  m 

*0  c 

X  ft 


C  3  »- 

X  w  w  ft  X  IS 

>  2  Zl  HH  • 
<  ft  S4  »-  ♦  2 
*-  ^  wh  ZZ  w 
lu  k  k  W  O  ■  k 
CHHSOZ W 


•  **22 


w  w 

a  s 

*s? 


-1  ©  ••  0- 
»  z  **  w»-  ft  tel 
ftia >  au4  js 
i-x 

J  »-  »-  IBM* 
klliiz  MffO 
ftauOik% % 
♦  ♦ 


>  -iff- 
ft  «0 


C  ft  9 

cuu 


"E 

-!■ 
ft  K 


::s 

»-  ~a 

00  tel  MO 
«»-► 
0*1  X 
I  or  O 
►  I  O 


o  e  o 

—  m  m 

u  u 

=>  C  O  O  C  CJ  o  o  o  o  c  c  ^coooccoooc  oo  ooo  oo  ooo  ooooco 
*  c  o  o  o  c  c  c  o  e  oo  c  coeeoeocoooooooeeocooeeooo 
c  ^ooococinoeecoooec  oocoooc  jinoeoooir<eKcaoo 

•.«.  I,  ^  \  v  ^  K  *.  k  K  • 

ir  tMOftc  MfuM  i\iki*  i/i  *  e~«At»A9ir4t<*c  oo  om 

9ccer(09(r999909ft999cococeooeeooo*-«**>««*«M«**«M«<«**« 

ir  ft  ►  r  ft  c  »  c  mmo*  tr  o  h  •  9  o-aw  o  m  r* 

00  «  C  10  9  9  9  0  O  9  9  0  9  9  COOOOOoOOO—  —  <n  *•  Ally  MV  M  M  ftt  IV 


182 


•  .000  HRITEOOO.OI) 

7.000  01  IM)*N*T(« 


X  M 

N. 

1*1  1*1 

K  C  - 

C 

«>  nz  hm 

4  ft  4. 

• 

J  Jff  X  X 

X  ft  *V 

44  3a  i 

K  *H 

•x. 

ft  O.K  X  X 

lb  wr» 

X 

•4 

•  •4404  •  *l*|  44  44 

I»i»iO 

ft  ft 

ItIUZX  ■  ■  ft  ^  ^ 

•X 

U  H  «  « 

2  ft 

C 

0  0  ft  ft  — ft»  ^  (Vil4 

fti  ft 

31*4 

•4  r 

*- 

J_4XVOClMww 

•5  V 

a©  ft  ~  ■ 

K  ►> 

Q 

4  4«4  r>2  ft  o  2  ft  « 

-J  —• 

U«K 

ft»  1*1 

e 

u 

ft  ft  -4fti~^ft  or  ft  U 

ft  ►* 

O  ft  uJ 

ft 

*  ft 

4- 

f  tKK  M  M  »X  >  K 

ft  ^ 

1*1  _!»-•* 

o 

%  #-* 

© 

■  I  ww44»-  «4^MK 

•  ft  1 

ft  U  U  4 

ft  ft 

V  ft 

• 

4 

—  IM2  2  X  31  C  K>—  X  4 

m  *4 

M  H(9 

JK  4 

O  * 

►- 

© 

C  C  ft  ft  Z  V  vft  0 

—  «  i 

k.  2  C  * 

*WW 

J>  o 

e 

*- 

2  2  00K  K  *2  2\K  ^ 

ft  ft 

UIC2» 

ft  X  o 

4> 

• 

^^44  W  0  X  ft  *•»  ft-  U 

ft  « 

6hm| 

S 

ft  4*-  2 

IV.  •  ru 

IV  X 

4 

C  C  «  «00  *ww#4  4*> 

44  m 

•-  X  C  U.  X  4  4  ^  ft 

•»«  a  h  i  \  2(v«r  j  k  • 

z  v  a  or  c  w*.  »orco 

-JO-/  •  C  2Ui\2wa  «kKhi 

W  ft  J  C  (9  U.  Jl  •2*>**‘OhJC« 

ft  *-•  »  *-_iftv"&«rvi5ikc* 

^  ■  i  •  ft.  ft  ►*  c  c  •*  w ,-  **  ** 

■  7 I H I 222N2COCC 


■  ft  X  ►* 

:  %  4i  ►*  \ 

■  vk  K  ft 

IK  U«Z  =  UiZ 

4  &  a  c  e  c 

*  r*  a  ft  ? 

X  *  K  Off  K  Z 

C  *  W  2  I?  z  c 

li  O  K  lU^H  W 
♦ 


N  z  c  o  ©  c 

*■**••••  ft) 

a  w  c  ft 

ft-no  JOtel  ?  7 
x  *  *  >  •  &  h 
ft  ft  -i  -j  > r  x  ft  ►* 


coco 
c  c  c  c 

©CPC 

.  4  • 

©  O  C  - 

—  —  rv  rv. 


c  r,  o  c  c 
o  ©  ©  c  c 
c  5  cr  / 

«.«.«•  i, 
•W  ft  4T  IT 
A»  \  i 


iUJ^CmmkhhIu 


Off  C  C  4 

»  e  r  c  e  i 

id©  w  —  « 

•.«.«.*.  t. 

iff  c  - 
l  K  P*  f*  W  I 


••  •  •  •  *0  Z  2«*4  % 
K  **n  ••  -OOKKZ  4  4  WK  «Q 
uiww  •  «C«ft  04  k  K1  u;k 
X  2  ZK  k  -  •  •  *44  ff  Z  * 

KSBCOffa^rs  «  «  VK  K 

*>'«  *  •  «4^  ».  •  *JP  V  W 

•  »Z«tt^^KKKhCCOMZ2 
Ui  MJ  D  K  K  *4  rv/fti  lb  u  UI2  <«•?: 

:  OUKKKKfff  CC  *«-w«ftft< 
ft  ft  O  *  «  w  v  AW  «.».  ff  ■  «  cu  I 
K  Kff  •  •■©  C Q«Mf>«93  «-«V  •«  48  ’ 

•  •  RIM  ZOOCff  21&.C  < 

HNJI  1447244  IK  Kff  J  .>1 
-/V-l  0  ■  WWW  ww  w  w  J  .J-  Jl 
CC4  4>Mtlklklb^kU.IKUl44  I 

i*  ZUkKmmmhvhmmCOKUOI 


>  c  O  ©  C  G  C  O  O  4 
;  o  o  c  c  c  eo  e  < 
;  o  o  o  c  c  0001 

i.  «.  i  4  •■  i  •.  «.  •. 
Jftl  4  N  tff  o  « 

>  tr  *\9t  **•*  tj  i 


:  >  c  o  o  o  c  < 
csccef < 

■  C  C-  o  o  a  e  ■ 

•.  b  <.  «.  «  «.  «. 

<•4  4IT  ff  KC  c 

I  ft  3  ft  ft  O  ft  « 


9 

c 

Z 


o 

* 


o 

z 


•3 

tel 


C 

2 


•  •  •  • 
w  mu  mi 
m  mm  m 


»«3 
l  -i  t* 
I  9  C 
.  k  Z 


■  9  m 

—  MM 

o  o  o 


•  k- 

9  O 
O  Z 


■  *> 
O  f 

z  • 

^  o 

O  flp 
•  •«  ■ 


mi  k- 

«vi  * 


•  D 
*•  O 

•  9 
H-  « 
O  • 


HOC  O 

•w  «c  z 

CD  —  O  • 
H  •« 

H  H  «  • 
•O  •  O 
9  •  O  Z 


kki 

(9  (JO 


•  & 
h  a 

C  H 
•  mi 
~  r 

U  H  k- 


CSC 

2 


m  a 
e  ♦ 

•  * 
a  ft 


z  z 

H  H 
■  ■ 

Al 
O  (9 


j  «m 
mi  im  a 

o  w  w 

«  z  z 
a  3  3 

A.  to.  m 

HU  u 

mi  *o  <c 


x  a. 
u  & 


C  W 

X  o 


ft  _i  -1 ' 

9  9  9  1 

(5U  U> 


Ulu  •"  • 

AHR 

^HM** 

*  Ni  O 

Z 

UHO  • 
U«^H 

•  «MC 

ac  z 

-J  9  2 


mu  (J 

a  i  ^  • 

a.  9  o  9 

<J  H-  <3  H- 

tu  *  taJ 

'  •!  C^X 
an  c  ch 

O  «  «tf  4t 


•  o  Z  Z  wft 
»  *03X3 


C5 

f 


U 

<3 


(3  CH 

•  u  ♦ 
m  ♦  • 
tt  civ 

W  H  V 

Z  ~  Z 
3X3 


9  U 
H  (9 

mi  % 
X  4D 
H  19 
*  * 
•  9 

me 


4»* 

% 

%  w 


m  mm  m 

z  z  z  z  • 

»  X  • 

•  H 

•H  Z 

•  •  •  • 

n  r« 

%  HI 

mi  z 

mix  k 

«  ■  «  ■ 

C  O  O  O  • 

♦  9 

Oft 

Cl  3 

**mc  9 

•  •  •  «o 

♦  H%  Z 

H 

-1  x 

J  X  H- 

oooo 

M  M«4  »^2 

*  mi  •  O 

X 

9  *"» 

9  r%UJ 

z  z  z  z  ^ 

•  •  •  «9 

♦  XX#-. 

•  ft 

m  o 

k  u  cr 

rs  ^  ^  (j 

4J 

U 

U 

HHH  H  • 

«►  »k 

mi  x 

•  k- 

•H  « 

«mc  9  H- 

k- 

H* 

1- 

oooo 

♦  »o  •-* 

•j 

«  ft  *- 

•  4-  (V 

o  oo  u  * 

* 

* 

% 

•  •  •  •  m 

«  ^  #c  « 

— i  ►- 

u  ^  x 

mso 

•  •  •  #c 

*> 

m 

40 

—  M9I  9  O 

♦  •*»  *o 

9  *** 

H  ^  OU  « 

ou  z 

HHH  HH 

k- 

k- 

k- 

HHHHI 

*  4 

Ik  ? 

*U  Z  H  U  2  H  • 

J  J  J  J  • 

% 

* 

UUUU  • 

♦  •* 

CH«wH 

19 

M 

9 

♦  ft  ft  z 

a  x  x 
coo 


0) 

e  •  •  • 

♦  900 


WZ^MJ  JJ 


-«IJ  JJ 

o  o  •  •  • 

Z  Z  —  IM  fC 
•  ■  19  OO 
CRwww 

com.  mm. 

2  2  H  NH 


J  U  U  KU 

•  m  cue 


U  J  J  H  J 
vJJUJ 


U  HU 
C  >  C 
#V 

-I  ►»  J 

J  U  J 


h*  u  uk  U  UiL  U  I 


> -J  -J 

:  9  «* 

.  u  u 


luuu 

9  9  9  9 
Hk  k  k 

O  w  w  w 


.  C  « 

z  ♦ 
•  •  « 
X  -J  41 
O  *  *  ♦ 

•  o  m  4i 

z  *  « 

•  9  r  41 

O  •  ♦ 

•  om 

H  p«- 

J  c  «% 

•  z  <e  x 

9  •  C  X 

HH  -  w 
UO  w»- 
9  Z  U/  9 
k  •  4-  X 

WVH  1 

m  m  a  c 

H  H  S  k 


~  »  9 
2  -J 

>  Jk 
9 

C  X 

mo 

%z 


HI  Z 
m  a  x 

ft  c  o 
k  u.  m 


C  CO 
—  fM  — 


o  o  c  o 
o  o  c  c 
e  c  9C 
.  1  I.  •.  I. 
9  tfl  ft  H> 
IT  IT  U*  LT 


ooocoooooc  ceooooco  oo  o  c  o  c  c  oooocoe  ooc  oooo 
o  o  e  o  c  c  o  cc  c  ec-mcoif  c  oooo  ccoccoecocoor  ooco 
ooooecoooc  o  —  —  — •  ~  —  c.  oooo  ez  o  c  o  ooo«mc4iroococ 

k  X  «.  k  •  k  •.•«*.  k  «.  I.  k  •.  •.  k  t.  «.  ««.•.«  I.  «.!«.(  k  •.  k  • 

o**m  cq  to  ■©h-ecooooo  o  cornea  it 
IT  11)44  44  ftftftftftftftftft  HH  H  H««t  «  O  «  ft  9  (XI  « 


v< 

Al 

X 

X 

u 

o 

ft 

ft 

4 

4 

ft 

k. 

44 

4\ 

% 

4 

— 

A 

44 

4* 

& 

ft 

— 

Vi 

a 

a. 

ft 

XI 

& 

a 

o 

c 

ft 

"N. 

ft 

a 

w 

W 

4* 

• 

u 

e 

44 

44 

mm 

w 

2 

2 

Z 

44 

4 

* 

4 

• 

2 

X 

— 

nj 

XI 

2 

V  4 

4 

a. 

ft 

ft 

*- 

•  — 

<v 

& 

ft 

4 

a 

*4  a 

a 

a 

c 

— 

«->  & 

a 

4 

4 

ft 

•  19 

13 

44 

44 

ft 

% 

2  4 

4 

— 

A 

w 

4 

44 

44 

X 

X* 

M 

ft 

a  — 

Al 

to 

to 

z 

x 

X 

IS 

ft 

4 

« 

4  r>  4) 

bi 

• 

1 

ft 

% 

•  c»  ft 

® 

w 

w 

ft 

w 

ru  4  ~ 

w 

44 

44 

4 

H  4  M 

— 

X 

z 

— 

— 

44 

44  % 

•  —  z 

2 

4 

4 

— 

44  — . 

— 

•  4 

4  >-  4 

4 

— 

A* 

— 

X>  XI 

ft 

u 

w  K  — 

A 

X 

4- 

4 

«.  4 

44 

►4  4- 

4- 

to 

to 

44  — 

—  44  x> 

— 

-1  * 

4  Wi  14 i 

to 

& 

* 

—  W 

M  4»  XI  44 

«  X 

4  — 

»z  a 

® 

■ 

• 

4  ft 

z  or  4ft 

X  4 

2  to 

13  to  4 

4 

4 

4 

Al  * 

a.  *  Ai  « 

ft 

8  U 

4  *K  « 

44 

— 

44 

—  44 

^44  W  44 

X  ft 

O  2 

V  4  2  — 

Al 

K\ 

ft 

ft  XI 

ft  AJ 

44  4 

2  |L> 

4-N  C3 

13 

13 

19 

44 

44  «  W 

X  4-  «  W 

X-  44 

C  •  IL  ^ 

w 

44 

*r 

XI  44  OQ 

^  ft  » 

X  XI 

m  h 

«  J  44 

44 

4* 

.  4- 

• 

w  K'  ft 

4-  ft  —  ft 

—  4 

*-«  o 

•  c3  ♦  a 

XI  A 

z 

z 

♦ 

•5  44W  4- 

X  -X  W  ♦ 

X  4. 

x  2 

2  2  4  ft  4 

•  4 

mm  m> 

XI  • 

4* 

ft  X  «  4 

44—  ft  — 

W  p 

ft  x* 

4  (3  4  — 

4  a 

4X1 

4  ft 

ms 

4  4  ft  — 

2  XI  ft  XI 

2  •-» 

*  %*«(» 

a*o 

*A19 

ft  15 

— 

mm 

—  44  a  4 

ft  4  ♦  4 

ft  4 

Al  14 

4-4  *4  (4 

44  « 

4*  4% 

<■4  ^ 

X 

X 

XI  X  **  ^  XI 

X»  XI  44  fU 

X  — 

C 

•  x.  -•  e  © 

o  © 

©  © 

©  © 

w 

mm 

*9  w  Al  4- 

ft  44  —  XI  W 

ft  — 

* 

4X  ft 

©  *4 

o  4 

o  — 

44 

X- 

XI  tS  Wi  4  ft 

•  ft  X  4ft 

*  44 

—  — 

mm  mm 

«4  — 

2 

2 

w  £  ay  xi  « 

•  ft 

ft  -J 

X«  P  X  4  % 

4  4 

4  4 

4  4 

4« 

ft  1*  4W  — 

©  —  X  W  44 

©  IU 

o  o 

4UUCC 

®  « 

•  ft 

ft  ft 

2 

• 

ft 

4ft  *4  ft  — 

«  —  —  ft  XI 

®  4 

M  2 

*1^00 

©  o 

©  ©  XI  o  ©  «* 

O 

— 

© 

44  |f1  M  £  W 

—  W  »-  «  4- 

—  — 

X  4* 

4*  ft  — 

V*  — 

mm  mm 

X 

«  4  K 

X- 

X 

ft 

ft 

XI  — 

«  ft  t  44ft 

P  4- 

« 

ft  %»4# 

w  w 

W  44  (J  V  (J 

X* 

44 

ft 

ft  ft  w  «  ft  —ft 

—  ft  w  fV  ft 

A  4» 

to  « 

«  a  >iu  ft 

to  to 

to  to 

XI 

to  to  ft 

U 

* 

ft 

V 

X  ft  —  ft  ^  I 

ft  •  2  w  1 

ft  © 

> 

4  KH  K 

X  X 

►  X* 

a 

Xh  ft 

2 

© 

ft 

X-  ft  *»ft  w  0>  44 

44  44  *1  ft  44 

44  ft 

(9  (9  M  H 

— 

44  44 

A 

44  H 

■« 

H  44  « 

3 

• 

ft  ft 

C  XI  42  ft  — 

•XI  ft  ft  ft 

•  4 

%  n  ft  ft  h  ft  a:  4- 

«  <fc 

ft  ft  a:  ft 

ft 

—  #4 

c 

2 

*fft2X4«3«  * 

O  4  ft  |  4 

o  m 

M  4 

•  «2  t  KUT  f  USB 

y  it  w 

X  X 

* © 

•<• 

X 

x  o  *■•  ■  ft  ^  <k 

.  •  Ai  *  44  — 

•  © 

X  w 

ft 

mm  mm 

ft 

<■»  4-»  ft 

ft 

*  w 

w 

X* 

H  I  H  u  IV  » 

ft  W  •  #4  w 

a  — 

ft  X 

4  JK  -  - 

4  a  «  *  * 

ft 

X«X2 

3 

2 

zm  a  «  4  a: 

■ 

2ft©  4ft 

« 

2  44 

41 

4-  a  2  o  o 

«  o  c  « 

30  XI 

C  O  ft  ft 

(A 

£ 

2 

ft 

© 

O 

C  2  ©  <bvo 

•  ■  «  XI  W 

© 

•Iftl 

•  2 

UJ  M  2  ?  4 

2  ?k 

2  2 

x  z  2  x  3 

ft  -3ft 

4*  o 

ft 

2 

a  ft  ft  —  -j  w  i 

9XX A  W  £ 

»xx 

(©  L9  ww|U 

w  w 

ft 

W  V(J 

vwOXOXX  ■ 

■ 

xoc 

2 

fix  ft  ft  • 

• 

W  P  ft  M 

■ 

44  M 

M  C 

4  8  IL 

a 

14  14 

a 

a  a 

X 

ft  ft  ft  ix  a 

r  2 

ft  4- 

44 

ft  j 

13  0 

o»-»  ft  oft  m  — 

— 

ft  AI  ■  — 

XI 

ft  ft 

►*  kft. 

4 

—  — 

• 

44  X 

ft 

44  44  a  ft  n 

U  X4 

ft  X 

2 

ft  b 

J  •> 

U 

UOQOUXX9 

—  ft  X  *- 

ft 

x*  C 

♦ 

♦  ♦  ♦ 

♦ 

♦ 

♦ 

♦ 

©  e  OOOOOCOOOOOOC  OOOOOJCOOOOCCO  OOOCOGQOOC  COO© 

ooccecoooc  ececooecoccircoocooo  oeccc  c  o  oc  c  c  o  o  c 
or>o»roooirociroc.iro  ©ir©©-*ftAiiAarfi©©o©  oocoocoec  c  aoc  o 

.  «.  •.  I.  i  «,  «  f,  «  «.  «.  «.  I,  ^  i  «L  •.  1  I.  ^  t,  t.  «  l,  ^  4  • 

£  h-  <  99  9  • 

c  «  4  «  c  £  a  a  ©  o  a  a  a  a  ©  ©  »  a  a  a  a  o  a  a  a  ©  c-  o  o  cocooe*-*<*-*-  _  —  M  — 


a  a.  a  a  f\»a  a  m\  a  a  a  <v  a  n.  a  ft  ft 


ELECTRONICS  TECHNOLOGY  AND  DEVICES  LABORATORY 
MANDATORY  CONTRACT  DISTRIBUTION  LIST 


101  Defense  Technical  Information  Center 
ATTN:  DTIC-TCA 
Cameron  Station  (Bldg  5) 

012  Alexandria,  VA  22314 

203  GIDEP  Engineering  &  Support  Dept 
TE  Section 
PO  Box  398 

001  NORCO,  CA  91760 

205  Director 

Naval  Research  Laboratory 
ATTN:  CODE  2627 
001  Washington,  DC  20375 

301  Rome  Air  Development  Center 
ATTN:  Documents  Library  (TILD) 

001  Griffiss  AFB,  NY  13441 

437  Deputy  for  Science  &  Technology 
Office,  Asst  Sec  Army  (R&D) 

001  Washington,  DC  20310 

438  HQDA  (DAMA-ARZ-D/Dr.  F.D.  Verderame) 
001  Washington,  DC  20310 

482  Director 

US  Army  Materiel  Systems  Analysis  Actv 
ATTN:  DRXSY-MP 

001  Aberdeen  Proving  Ground,  MD  21005 

563  Commander,  DARCOM 
ATTN:  DRCDE 

5001  Eisenhower  Avenue 
001  Alexandria,  VA  22333 

564  Cdr,  US  Army  Signals  Warfare  Lab 
ATTN:  DELSW-OS 

Vint  Hill  Farms  Station 
001  Warrenton,  VA  22186 

705  Advisory  Group  on  Electron  Devices 
201  Varick  Street,  9th  Floor 
002  New  York,  NY  10014 

579  Cdr,  PM  Concept  Analysis  Centers 
ATTN:  DRCPM-CAC 
Arlington  Hall  Station 


001  Arlington,  VA  22212 

602  Cdr,  Night  Vision  &  Electro-Optics 
ERADCOM 

ATTN:  DELNV-D 
001  Fort  Belvoir,  VA  22060 

603  Cdr,  Atmospheric  Sciences  Lab 
ERADCOM 

ATTN:  DELAS-SY-S 

001  White  Sands  Missile  Range,  NM  88002 

607  Cdr,  Harry  Diamond  Laboratories 
ATTN:  DELHD-CO,  TD  (In  Turn) 

2800  Powder  Mill  Road 
001  Adelphi,  MD  20783 

609  Cdr,  ERADCOM 

ATTN:  DRDEL-CG,  CD,  CS  (In  Turn) 

2800  Powder  Mill  Road 
001  Adelphi,  MD  20783 

612  Cdr,  ERADCOM 
ATTN:  DRDEL-CT 
2800  Powder  Mill  Road 
001  Adelphi,  MD  20783 

680  Commander 

US  Army  Electronics  R&D  Command 
000  Fort  Monmouth,  NJ  07703 
1  DELET-MQ 
1  DELEW-D 
1  DELET-DD 

1  DELSD-L  (Tech  Library) 

2  DELSD-L-S  (STINFO) 

34  Originating  Office 

1  DELET-MF 

681  Commander 

US  Army  Communications  R&D  Command 
ATTN:  USMC-LNO 
001  Fort  Monmouth,  NJ  07703 
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103 

Code  R123,  Tech  Library 

475 

Cdr,  Harry  Diamond  Laboratories 

DCA  Defense  Comm  Engrg  Ctr 

ATTN:  Library 

1800  Wiehle  Ave 

2800  Powder  Mill  Road 

001 

Reston,  VA  22090 

001 

Adelphi,  MD  20783 

104 

Defense  Communications  Agency 

477 

Director 

Technical  Library  Center 

US  Army  Ballistic  Research  Labs 

Code  205  (P.  A.  Tolovi) 

ATTN:  DRXBR-LB 

001 

Washington,  DC  20305 

001 

Aberdeen  Proving  Ground,  MD  21005 

206 

Commander 

*481 

Hairy  Diamond  Laboratories 

Naval  Electronics  Laboratory  Center 

ATTN:  DELHD-RCB  (Dr.  J.  Nemarich) 

ATTN.  Library 

2800  Powder  Mill  road 

001 

San  Diego,  CA  92152 

001 

Adelphi,  MD  20783 

207 

Cdr,  Naval  Surface  Weapons  Center 

482 

Director 

White  Oak  Laboratory 

US  Army  Materiel  Systems  Analysis  Actv 

ATTN:  Library  Code  WX-21 

ATTN:  DRXSY-T,  MP  (In  Turn) 

001 

Silver  Spring,  MD  20910 

001 

Aberdeen  Proving  Ground,  MD  21005 

314 

Hq,  Air  Force  Systems  Command 

507 

Cdr,  AVRADCOM 

ATTN:  DLCA 

ATTN:  DRSAV-E 

Andrews  Air  Force  Base 

PO  Box  209 

001 

Washington,  DC  20331 

001 

St.  Louis,  MO  63166 

403 

Cdr,  MICOM 

511 

Commander,  Picatinny  Arsenal 

Redstone  Scientific  Info  Center 

ATTN:  SARPA-FR-5,  -ND-A< 

ATTN.  Chief,  Document  Section 

-TS-S  (In  Turn) 

001 

Redstone  Arsenal,  AL  35809 

001 

Dover,  NJ  07801 

406 

Commandant 

515 

Project  Manager,  REMBASS 

US  Army  Aviation  Center 

ATTN:  DRCPM-FFR-TM 

ATTN:  ATZQ-D-MA 

001 

Fort  Monmouth,  NJ  07703 

001 

Fort  Rucker,  AL  36362 

517 

Commander 

407 

Director,  Ballistic  Missile  Defense 

US  Army  Satellite  Communications  Agcy 

Advanced  Technology  Center 

ATTN:  DRCPM-SC-3 

ATTN:  ATC-R,  PO  Box  1500 

001 

Fort  Monmouth,  NJ  07703 

001 

Huntsville,  AL  35807 

518 

TRI-TAC  Office 

418 

Commander 

ATTN:  TT-SE 

HQ,  Fort  Huachuca 

001 

Fort  Monmouth,  NJ  07703 

ATTN:  Technical  Reference  Div 

001 

Fort  Huachuca,  AZ  85613 

‘For  Millimeter  &  Microwave  Devices  Only 
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519 

Cdr,  US  Army  Avionics  Lab 

608 

Commander 

AVRADCOM 

ARRADCOM 

ATTN:  DAVAA-D 

DRDAR-TSB-S 

001 

Fort  Monmouth,  NJ  07703 

001 

Aberdeen  Proving  Ground,  MD  21005 

520 

Project  Manager,  FIREF1NDER 

614 

Cdr,  ERADCOM 

ATTN:  DRCPM-FF 

ATTN.  DRDEL-LL,  -SB,  -AP  (In  Turn) 

001 

Fort  Monmouth,  NJ  07703 

2800  Powder  Mill  Road 

521 

Commander 

001 

Adelphi,  MD  27083 

Project  Manager,  SOTAS 

617 

Cdr,  ERADCOM 

ATTN:  DRCPM-STA 

ATTN:  DRDEL-AQ 

001 

Fort  Monmouth,  NJ  07703 

2800  Powder  Mill  Road 

531 

Cdr,  US  Army  Research  Office 

001 

Adelphi,  MD  20783 

ATTN:  DRXRO-PH  (Dr.  Lontz) 

619 

Cdr,  ERADCOM 

DRXRO-IP  (In  Turn) 

ATTN.  DRDEL-PA,  -ILS,  -ED  (In  Turn) 

PO  Box  12211 

2800  Powder  Mill  Road 

001 

Research  Triangle  Park,  NC  27709 

001 

Adelphi,  MD  20783 

556 

HQ.  TCATA 

701 

MTI  —  Lincoln  Laboratory 

Technical  Information  Center 

ATTN:  Library  (RM  A-082) 

ATTN:  Mrs.  Ruth  Reynolds 

PO  Box  73 

001 

Fort  Hood,  TX  76544 

002 

Lexington,  MA  02173 

568 

Commander 

703 

NASA  Scientific  &  Tech  Info  Facility 

US  Army  Mobility  Eqp  Res  &  Dev  Cmd 

Baltimore/Washington  Inti  Airport 

ATTN:  DRDME-R 

001 

PO  Box  8757,  MD  21240 

001 

Fort  Belvoir,  VA  22060 

704 

National  Bureau  o'  Standards 

6C4 

Chief 

Bldg  225,  RM  A-331 

Ofc  of  Missile  Electronic  Warfare 

ATTN:  Mr.  Leedy 

Electronic  Warfare  Lab,  ERADCOM 

001 

Washington,  DC  20231 

001 

White  Sands  Missile  Range,  NM  88002 

.  707 

TACTEC 

606 

Chief 

Batelle  Memorial  Institute 

Intel  Materiel  Dev  &  Support  Ofc 

505  King  Avenue 

Electronic  Warfare  Lab,  ERADCOM 

001 

Columbus,  OH  43201 

001 

Fort  Meade,  MD  20755 

AD^AUO  6*3  MOTOROLA  INC  SCOTT  JO  ALE  AZ  60VERNMENT  ELECTRONICS  DIV  f/% 

DOUBLY  ROTATED  CUT  SAM  DEVICES. (U) 

JAN  62  D  F  WILLIAMS.  F  Y  CHO  DAAK20-70-C-027S 

UNCLASSIFIED  DELET-TR-T9-02TS-F  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BIIRLAU  Of  SIANOAHPS  1%.<  A 
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Coordinated  Science  Laboratory 
University  of  Illinois 
Urbana,  Illinois  61801 

ATTN:  Or.  Bill  J.  Hunsinger  (1) 

Dr.  J.S.  Bryant 
OCRD 

ATTN.  OARO-ARP 

Washington,  DC  20310  (1) 

Dr.  R.  LaRosa 

Hazeltine  Corporation 

Greenlawn,  New  York  11740  (1) 

General  Electric  Co. 

Electronics  Lab 
Electronics  Park 
Syracuse,  NY  13201 

ATTN:  Mr.  S.  Wanuga  (1) 

Air  Force  Cambridge  Labs 
ATTN.  CRDR  (Dr.  P.  Carr  &  Dr.  A.J. 
Slobodnik) 

Bedford,  MA  01730  (2) 

Mr.  R.  Weglein 

Hughes  Research  Laboratories 

3011  Malibu  Canyon  Road 

Malibu,  California  90265  (1) 

Mr.  H.  Bush  CORC 
RADC 

Griffiss  Air  Force  Base 

New  York  13440  (1) 

Or.  Tom  Bristol 
Hughes  Aircraft  Company 
Ground  Systems  Group 
Bldg  600/MS  0235 
1901  W.  Malvern 

Fullerton.  CA  92634  (2) 

Commander,  AFAL 

ATTN:  Mr.  W.J.  Edwards,  TEA 

Wright-Patterson  AFB,  Ohio  45433  (1) 


Anderson  Laboratories,  Inc. 

1280  Blue  Hills  Ave 

ATTN:  Dr.  A.A.  Comparini 

Bloomfield,  Conn.  06002  (1) 

Mr.  Henry  Friedman 
RADC/OCTE 

Griffiss  AFB,  NY  13440  (1) 

Autonetics,  Division  of  North 
American  Rockwell 
P.0.  Box  4173 
3370  Miraloma  Avenue 
Anaheim,  CA  92802 

ATTN:  Dr.  G.R.  Pulliam  (1) 

General  Dynamics,  Electronics  Division 
P.O.  Box  81127 
San  Diego,  CA  92138 

ATTN:  Mr.  R.  Badewitz  (1) 

Texas  Instruments,  Inc. 

P.O.  Box  5936 

13500  N.  Central  Expressway 
Dallas,  Texas  75222 

ATTN:  Dr.  L.T.  Clairborne  (2) 

Raytheon  Company 

Research  Division 

28  Seyon  Street 

Waltham,  Massachusetts  02154 

ATTN:  Dr.  M.B.  Schulz  (1) 

Sperry  Rand  Research  Center 
100  North  Road 

Sudbury,  Massachusetts  01776 
ATTN:  Dr.  H.  Van  Oe  Vaart  (1) 

Microwave  Laboratory 
W.W.  Hansen  Laboratories  of  Physics 
Stanford  University 
Stanford,  CA  94305 

ATTN:  Dr.  H.J.  Shaw  (2) 
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Polytechnic  Institute  of  Brooklyn 

Advanced  Technology  Center,  Inc. 

Route  No.  110 

Subsidiary  of  LTV  Aerospace  Corp. 

Farmingdale,  NY  11735 

P.O.  Box  6144 

ATTN:  Dr.  A.A.  Oliner 

(D 

Dallas,  Texas  75222 

ATTN:  Mr.  A.E.  Sobey 

(1) 

Westinghouse  Electric  Corp. 

Research  &  Development  Center 

United  Aircraft  Research  Labs 

Beulah  Road 

ATTN:  Dr.  Thomas  W.  Grudkowski 

Pittsburgh,  PA  15235 

East  Hartford,  Conn.  06108 

(D 

ATTN:  Dr.  J.  DeKlerk 

(1) 

Science  Center 

Stanford  Research  Institute 

Rockwell  International 

Menlo  Park,  CA  94025 

Thousand  Oaks,  CA  91360 

ATTN:  Dr.  A.  Bahr 

(1) 

ATTN:  Dr.  T.C.  Urn 

(D 

International  Business  Machines  Corp. 

University  of  Southern  CA 

Research  Division 

Electronic  Science  Lab 

P.O.  Box  218 

School  of  Engineering 

Yorktown  Heights,  NY  10598 

University  Park,  Los  Angeles 

ATTN:  Dr.  F.  Bill 

(1) 

California  900 

TRW 

ATTN:  Dr.  K.  Lakin,  SSC  303 

(1) 

Defense  and  Space  Sys  Group 

SAWTEK,  Inc. 

One  Space  Park 

P.O.  Box  7756 

Redondo  Beach,  CA  90278 

2451  Shader  Road 

ATTN.  Dr.  R.S.  Kagiwada 

(D 

Orlando,  Florida  32854 

ATTN:  Mr.  S.  Miller 

(1) 

Tektronix  Inc. 

P.O.  Box  500 

Prof.  P.C.Y.  Lee 

Beaverton,  OR  97077 

School  of  Engineering 

ATTN:  Dr.  R.  Li 

(1) 

Princeton  University 

Dr.  Fred  S.  Hickernell 

Princeton,  NJ  08540 

(D 

Integrated  Circuit  Facility 

Mr.  John  A.  Kusters 

Motorola  Government  Electronics  Oivision 

Hewlett-Packard 

8201  East  McDowell  Road 

5301  Stevens  Creek  Boulevard 

Scottsdale,  AZ  85257 

(1) 

Santa  Clara,  CA  95050 

(1) 

Prof.  H.F.  Tiersten 

Dr.  Tom  Young 

Jonsson  Engineering  Center 

Sandia  Laboratories 

Rensselaer  Polytechnic  Institute 

P.O.  Box  5800 

Troy,  NY  12181 

(D 

Albuquerque,  NM  87185 

(D 

McGill  University 

Dr.  William  J.  Tanski 

ATTN.  G.W.  Farnell 

Sperry  Research  Center 

Montreal  110,  Canada 

(D 

100  North  Road 

Sudbury,  MA  01776 

(1) 
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Dr.  B.A.  Auld 

Dr.  Robert  L.  Rosenberg 

E.L.  Ginzton  Laboratory 

Bell  Laboratories 

Stanford  University 

Crawfords  Corner  Road 

Stanford,  CA  94305 

(1) 

Holmdel,  NJ  07733 

(D 

Mr.  Marvin  E.  Frerking 

Dr.  B.K.  Sinha 

MS  137-138 

Schiumberger-Doll  Research  Center 

Collins  Radio  Company 

P.0.  Box  307 

855  35th  Street,  NE 

Ridgefield,  CT  06877 

(D 

Cedar  Rapids,  IA  52406 

(D 

Dr.  Robert  S.  Wagers 

Dr.  William  R.  Shreve 

Texas  Instruments,  Inc. 

HP  Laboratories 

13500  N.  Central  Expwy. 

1501  Page  Mill  Road 

P.0.  Box  225936,  MS  134 

Palo  Alto.  CA  94304 

(D 

Dallas.  TX  75265 

(D 

Dr.  Thomas  M.  Reeder 

Dr.  Richard  C.  Williamson 

MS  50-362 

Lincoln  Laboratory 

Tektronix,  Inc. 

P.O.  Box  73 

P.0.  Box  500 

Beaverton,  OR  97077 

(D 

Lexington,  MA  02173 

(D 

SUPPLEMENT  TO  DISTRIBUTION  LIST 

D.  Chrissotimos,  Code  763 

Army  Materials  and  Mechanics  Research 

National  Aeronautics  and  Space 

Center  (AMMRC) 

Administration 

Watertown,  MA  02172 

Goddard  Space  Flight  Center 

ATTN:  DMXMR-EO 

(D 

Greenbelt,  MD  20771 

(D 

Commander,  Picatinny  Arsenal 

Naval  Research  Laboratories 

ATTN:  SARPA-FR-S 

Code  5237 

Bldg.  350 

Washington,  DC  20375 

Dover,  NJ  07801 

(2) 

ATTN:  Dr.  D.  Webb 

(1) 

A.  Kahan 

HQ  ESD  (DRI) 

RADC/ESE 

L.G.  Hanscom  AFB 

Hanscom  AFB 

Bedford,  MA  01731 

(D 

Bedford,  MA  01731 

(D 

Commander 

Dr.  Robert  O'Connell 

US  Army  Missile  Command 

Department  of  EE 

ATTN:  DRSMI-RE  (Mr.  Pittman) 

Univ.  of  Missouri  -  Columbia 

Redstone  Arsenal,  AL  35809 

(D 

Columbia,  MO  65201 

(D 

Prot.  John  F.  Vetelino 
Dept  of  EE 
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Dr.  Robert  O’Connell 

Department  of  EE 

University  of  Missouri  -  Columbia 

Columbia.  MO  65201  (1) 

Prof.  John  F.  Vetelino 
Dept  of  EE 

University  of  Maine  •  Orono 

Orono,  ME  04469  (1) 

Dr.  Franz  Sauerland 
Transat  Corp. 

3713  Lee  Road 

Shaker  Heights,  OH  44120  (1) 

Mr.  Leon  Marchand 
Damon  Corp. 

80  Wilson  Way 

Westwood,  MA  02090  (1) 


